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ABSTRACT 


A  major  phase  of  this  report  involved  a  compilation  of  information 
and  critical  evaluation  of  that  information  as  it  pertains  to  the  diagnostic 
problem.  The  information  is  gathered  together  in  a  bibliography  of 
books,  reports,  and  published  articles  and  is  discussed  in  Section  II, 

The  investigation  of  the  feasibility  of  applying  the  theory  of  band  models 
to  the  diagnostic  problem  of  determining  temperatures  and  concentra¬ 
tions  has  led  to  the  development  of  a  mathematical  inversion  technique 
based  on  low-resolut.'on  emission-absorption  spectroscopy.  Band 
model  theory  and  mucp.  of  its  literature  is  summarized  in  this  report. 
The  inversion  of  the  Cu.^tis -Godson  approximation  of  a  statistical  band 
model  with  Lorentz  lines  of  constant  width  is  presented  in  detail. 
Experiments  were  performed  to  verify  the  technique  but  are  inconclu¬ 
sive  at  this  time. 
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Intensive  efforts  have  been  expended  In  developing 
non-interference  diagnostic  methods  based  on  spectroscopy  for 
determining  direct  in  situ  quantitative  values  of  the  chemical 
and  thermodynamic  properties  of  radiating  gases. 

The  radiant  energy  emitted  by  a  hot  gas  is  the 
integral  of  contributions  from  all  the  spatial  regions  within 
the  volume  of  the  gas  and  all  the  spectral  intervals  over 
which  it  is  active.  Spatial  integration  is  usually  possible 
but  spectral  integration  often  is  not.  To  predict  accurately 
the  radiation  from  hot  gases,  the  spectral  absorption 
coefficient,  k^,  as  a  function  of  frequency  must  be  known. 

And  so  despite  the  efforts,  there  still  exist  many  experimental 
situations  in  which  it  is  not  possible  to  obtain  sufficiently 
well  resolved  spectra  to  determine  the  information  desired. 
Furthermore  even  in  those  situations  in  whiph  the  resolution  is 
sufficient  to  define  the  individual  lines  of  a  vibration- 
rotation  band  k^  is  a  rapidly  varying  function  of  frequency  and 
the  computational  task  may  beccmie  quite  formidable  and 
expensive  when  integrating  over  frequency  even  for  a  homogeneous 
gas  of  simple  molecules  and  consequent  simple  spectrum.  Both 
for  interpreting  measurements  and  computing  numerical  problems, 
it  has  therefore  proved  desirable  to  develop  techniques  for 
averaging  over  fine  structure  in  the  spectrum.  The  averaging 
process  is  intended  to  smooth  out  the  lines  without  distorting 
the  band  contour  and  is  thus  referred  to  as  a  "band  model". 

Band  models  provide  the  averaging  processes  by  substituting 
for  the  true  spectrum  an  infinite  array  of  spectral  lines  of 
uniform  statistical  properties.  Four  types  of  band  models 
are  particularly  useful  and  will  be  briefly  discussed.  These 
are  the  Elsasser  model,  the  statistical  model,  the  random 
Elsasser  model,  and  the  Quasi-random  model. 

These  band  models  were  designed  to  allow  the  pre- 
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diction  of  the  infrared  spectral  intensity  of  an  homogensous 
gas  when  temperatures  and  pressures  were  known.  The  extent 
to  which  the  averaging  process  could  be  applied  to  infrared 
measurements  of  inhomogeneous  gases  with  gradients  of 
temperature  pressures  and  concentrations  required  investigation. 
The  Curtis-Godson  approximation  provides  a  technique  for 
calculating  the  radiance  and  transmittance  of  hot  inhomogeneous 
gases  by  determining  the  necessary  spectral  parameters  solely 
from  homogeneous  gas  data.  In  the  Curtis-Godson  approximation 
the  transmittance  of  the  inhomogeneous  path  is  obtained  by 
substituting  a  hypothetical  homogeneous  path  with  the  same 
transmittance.  The  homogeneous  path  can  then  be  treated  by 
established  band  model  methods. 

The  objective  of  the  present  program  was  fourfold: 

1.  to  acquire  an  understanding  of  band  model 
theory  and  to  present  the  information  and 
and  development  of  band  model  theory; 

2.  to  study  the  applications  of  band  model 
theory  to  various  physical  situations; 

3 .  to  examine  the  limitations  and  validity 
of  band  models  as  an  averaging  process 
for  determining  the  emissivity  and 
absorptivity  of  finite  spectral  bands 
containing  many  unresolved  lines  of 
water  vapor  and  carbon  dioxide  in  a 
radiating  inhomogeneous  gas;  and 

4.  to  determine  the  extent  to  which 

the  accumulated  knowledge  of  experimental 
techniques  and  spectral  parameters  of 
previous  investigators  could  be  used 
to  develop  a  diagnostic  method.  That  is 
a  method  based  on  band  model  theory  for 
evaluating  the  temperatures  and  partial 
pressures  of  hot  inhomogeneous  gases  from 
measured  spectral  intensities. 
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The  present  report  is  intended  to  be  first  of  all 
an  introductory  handbook  of  band  model  theory  and  second 
a  description  of  a  band  model  diagnostic  technique.  Much 
of  the  pertinent  literature  was  searched  and  reviewed. 
Information  was  compiled  and  evaluated  and  is  presented  as 
a  guide  to  the  literature  in  Section  II.  An  exposition  of 
the  theory  of  band  models  beginning  with  the  single  spectral 
line  and  extending  the  development  through  the  Curtis-Godson 
approximation  to  inhomogeneous  gases  is  set  forth  in  Section 
III.  .The  extensive  literature  search  led  to  a  restricting 
and  defining  of  the  diagnostic  problem  and  subsequently  to 
the  development  of  a  mathematical  inversion  technique  based 
on  low  resolution  emission-absorption  spectroscopy  for 
determining  temperatures  and  partial  pressures  of  infrared 
active  species.  This  is  treated  in  Section  IV.  Section  V 
describes  the  experiment  designed  to  test  the  consistency 
and  validity  of  the  computational  method. 

The  objective  of  the  experiment  was  to  seek 
spectral  emission  and  absorption  data  at  2.7|j,  in  the  radiation 
spectrvim  of  hot  water  vapor  produced  in  flames.  These 
experiments  to  provide  the  required  line-of-sight  radiance  and 
transmittance  data  were  carried  out  by  UTSI  personnel  in 
cooperation  with  the  ETF  Special  Projects  Group  utilizing  the 
special  burners  and  spectrometers  available  at  AEDC.  Section 
VI  contains  the  concluding  remarks. 
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SECTION  li 
,  LITERATURE  GUIDE 

A  major  phase  of  this  report  involved  a  compilation  of 
information  and  critical  evaluation  as  it  pertains  to  the 
diagnostic  problem.  This  information  is  gathered  together  in 
the  bibliography  and  will  be  commented  on  in  this  section  of 
the  report  and  referred  to  at  various  points  throughout  the 
report.  The  bibilography  is  a  collection  of  books,  reports, 
and  published  articles  which  have  been  gathered  together  for 
study  and  evaluation.  This  subject  is  weighted  down  by  a 
huge  unpublished  literature  in  the.  form  of  contract  reports, 
such  as  this,  which  do  not  become  readily  available  to  the 
general  scientific  community.  The  present  bibliography  con¬ 
tains  such  reports  because  they  sometimes  contain  information 
not  available  elsewhere,  but  the  emphasis  where  possible  is 
placed  on  books  and  articles  appearing  in  the  open  literature. 
The  present  collection  certainly  does  not  exhaust  the  subject 
even  in  the  English  language  and  hardly  touches  on  the 
foreign  literature  and  yet  it  is  believed  to  represent  much 
of  the  major  effort  in  this  area.  Although  many  approaches 
to  a  literature  commentary  suggest  themselves,  the  present 
commentary  will  for  the  most  part  be  directed  to  the  efforts 
of  groups  of  researchers  working  and  reporting  over  a  period 
of  several  years  because  that’s  how  much  of  the  work  has 
been  performed.  However,  the  accomplishments  of  individuals 
often  laid  the  groundwork. 

The  first  physically  satisfying  band  model  was  developed 
by  W.  M.  Elsasser,  in  1938  (Ref.  1)  and  in  1942  (Ref.  2)  to 
aid  in  predicting  the  transmission  of  Infrared  radiation 
through  the  atmosphere .  The  principle  of  the  random  or 
statistical  model  was  first  published  by  R.  M.  Goody  in  1952 
(Ref.  3).  These  two  models  provide  a  starting  point  for 
understanding  band  model  theory  and  remain  basic  to  the  study 
of  many  molecules . 

The  fundamental  physical  principles  of  band  model  theory 
are  most  completely  developed  in  the  very  excellent  book 
"Atmospheric  Radiation"  by  R.  M.  Goody  (Ref.  4).  This  book 
will  be  most  satisfying  to  those  who  have  some  background  in 
radiative  heat  transfer  and  spectroscopy  although  both  of 
these  areas  are  developed  from  basic  concepts  in  the  book. 

The  first  four  chapters  are  particularly  informative  and 
pertinent  to  radiation  theory  in  general  and  band  model  theory 
in  particular.  Each  chapter  is  capped  off  with  a  most  complete 
bibliography  wherein  many  references  are  accompanied  by 
remarks  suggesting  their  usefulness.  Many  excellent  diagrams 
and  graphs  enhance  the  exposition  in  the  text.  There  are 
fourteen  appendices  several  of  which  are  very  useful  tables. 

A  second  book  which  provides  considerable  background 
in  the  theory  and  techniques  of  infrared  radiation  and  the 
application  of  band  models  to  radiative  transfer  problems  is 
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the  Handbook  of  Military  Infrared  Technology  (Ref.  5).  As 
the  title  implies  the  book  is  exhaustive  in  coverage,  con¬ 
taining  hundreds  of  graphs  and  tables  and  references. 

A  very  useful  state-of-the-art  report  is  ’’Band-Model 
Methods  for  Computing  Atmospheric  Slant-Path  Molecular 
Absorption”  written  by  David  Anding  (Ref.  6).  This  report 
contains  a  summary  development  of  the  band  models  together 
with  many  tables  and  graphs  comparing  the  usefulness  of  the 
band  models  when  applied  to  the  constituents  of  the 
atmosphere.  The  report  by  Rochelle  (Ref.  7)  is  similar  in 
development  although  it  emphasises  the  application  to  thermal 
radiation  from  rocket  exhausts. 

The  theory  of  band  models  is  based  on  parameters  involved 
in  describing  individual  spectral  lines.  The  development  of 
spectral  line  theory  is  described  in  many  places  but-  the 
present  development  is  referred  to  the  two  papers  by  Benedict 
Herman,  Moore  and  Silverman  (Refs.  8  and  9)  and  the  appli¬ 
cation  by  Mason  and  Nielson  (Ref.  10).  In  extending  the  band 
model  theories  from  homogeneous  to  inhomogeneous  gases  the 
contribution  of  Curtis  (Ref.  11)  and  Godson  (Ref.  12)  are 
invaluable.  These  twelve  references,  together  with  the 
reported  work  of  G.  N.  Plass  provided  the  background  from  ' 
which  the  section  of  this  report  entitled  Band  Model  Theory 
was  written. 

One  need  not  look  far  into  the  literature  before  it 
becomes  evident  that  band  model  theory  was  developed  for  and 
applied  to  radiative  transfer  through  the  atmosphere.  A 
foremost  worker  in  this  area  has  been  Gilbert  N.  Plass  whose 
articles  have  spanned  nearly  twenty  years  (Refs.  13  through 
20) .  Much  of  his  effort  has  been  devoted  to  refining  band 
model  theory  and  examining  the  criteria  for  the  applicability 
and  limitations  of  the  different  models.  His  reported  work  is 
sufficiently  general  that  it  is  quite  as  pertinent  to  com¬ 
bustion  problems  as  to  atmospheric  transmission  problems. 

Plass  proposed  three  approximations  that  can  be  used  to  give 
relatively  simple  analytical  expressions  for  the  absorptance 
predicted  by  the  various  band  models  for . virtually  any  range 
of  pressure  and  concentration  of  absorbing  gas.  These  will 
be  discussed  in  detail  later  in  this  report. 

Plass  also  describes  laboratory  studies  appropriate  to 
the  selection  of  the  correct  band  and  the  correct  approxi¬ 
mation  for  the  particular  situation  of  Interest. 

A  very  active  experimental  group  orginating  at  Ohio 
State  consisting  of  D.  E.  Burch,  D.  Williams,  D.  Gryvnak,  and 
R.  R.  Patty  have  published  many  papers  on  the  results  of 
extensive  laboratory  measurements  of  the  total  absorptance  of 
the  major  near  infrared  absorption  bands  and  lines  of  gases 
present  in  planetary  atmospheres  (Refs.  21  through  33). 

These  gases,  in  particular  carbon  monoxide,  carbon  dioxide, 
and  water  vapor  are  often  the  gases  of  most  Interest  in 
combustion  phenomena  and  so  the  research  of  this  group  is  of 
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considerable  Interest  to  this  report.  In  these  experiments 
the  total  absorptance  vas  measured  as  functions  of  absorber 
concentration  and  total  pressure  providing  empirical 
relationships  for  a  considerable  range  of  these  variables. 

Of  particular  import  to  band  model  theory  is  that  Burch  and 
Williams  (Refs.  28)  examine  the  possibility  of  experimental 
realization  of  the  conditions  leading  to  the  Plass  approxi¬ 
mations  and  conclude  that;  (1)  laboratory  results  give 
satisfactory  approaches  to  the  weak  line  approximation  but 
not  to  the  strong  line  approximation;  and  (2)  it  may  often 
be  necessary  to  introduce  the  results  of  high  resolution 
studies  to  properly  select  appropriate  band  models. 

There  have  been  several  groups  active  in  the  application 
of  band  model  techniques  to  the  problem  of  radiative  transfer 
in  rocket  exhaust  gases.  Much  of  this  work  was  performed 
under  contract  to  the  National  Aeronautics  and  Space 
Administration.  A  summary  and  cross  section  of  that  work 
appears  in  a  publication  entitled  "Molecular  Radiation  and  Its 
Application  to  Diagnostic  Techniques"  edited  by  R.  Goulard, 
Purdue  University  (Ref .  34) .  This  is  the  record  of  the 
Proceedings  of  a  Specialists  Conference  held  at  the  Marshall 
Space  Flight  Center  Huntsville,  Alabama,  October  5-6,  1967, 

Two  research  groups  who  were  very  much  involved  in  the  NASA 
effort  and  whose  work  provides  the  beginning  point  for  the 
present  UTSI-AEDC  effort  were  the  General  Dynamics  Convair 
Group  and  the  Warner  and  Swasey  Group. 

The  General  Dynamics  Convair  Group  consisted  of  C.  C. 
Ferriso,  (deceased),  C.  B.  Ludwig,  J.  A.  L.  Thomson,  and  W. 
Malkmus.  Several  reports  and  journal  articles  (Refs.  35 
through  51)  document  their  research,  which  centered  on  the 
measurement  of  the  absorption  coefficients  and  fine  structure 
parameters  of  water  vapor  and  carbon  dioxide  at  temperatures 
between  1200°  and  3000°K.  These  data  were  then  incorporated 
into  a  practical  analytic  procedure  which  could  be  used  for 
the  evaluation  of  the  radiative  base  heating  of  rocket 
vehicles.  Their  experimental  setup  was  quite  unique  providing 
a  volume  of  high  temperature  water  vapor  above  the  combustion 
zone  of  a  long  three  inch  wide  flat  flame  burner.  Emission 
and  absorption  measurements  were  made  with  a  path  length  of 
twenty  feet.  Also,  nonisothermal  measurements  were  made  with 
four  five-foot  sections  at  different  temperatures  to  test  the 
nonisothermal  radiance  calculation  procedure  which  was  a 
modified  Curtis-Godson  method  worked  out , independently  by  the 
General  Dynamics  Convair  Group  and  the  Warner  and  Swasey  Group, 
This  calculation  was  the  basis  for  the  computer  program  used 
by  the  NASA  Marshall  Space  Flight  Center  to  calculate  the 
radiant  heating  of  a  rocket  base.  The  most  important 
consequence  of  the  General  Dynamics  Program  has  been  to 
generate  tables  of  water  vapor  absorption  coefficients 
covering  the  spectral  region  50  to  9300  cm"^  and  the  tempera¬ 
ture  range  300°  to  3000°K  and  tables  of  the  reciprocal  of  the 
average  line  spacing,  d,  covering  the  interval  1150  to 
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9300  cm“^  in  the  temperature  range  from  600°  to  3000°K. 

William  F.  Herget,  Rocketdyne  Division,  North  American 
Rockwell  Corporation,  whose  work  (Refs.  52  through  55)  is 
discussed  later  in  the  report,  contributed  smeared  high 
resolution  water  vapor  data  to  the  generation  of  the  above 
tables.  Fredrick  S.  Simmons,  Willow  Run  Laboratories,  The 
University  of  Michigan  has  also  been  involved  in  the  NASA 
effort  to  study  the  application  of  band  models  to  nonisothermal 
radiating  gases  (Refs.  56  through  59)  and  provided  comparison 
data  for  an  independent  check  to  the  General  Dynamics  data  for 
temperatures  greater  than  1000°K. 

The  Warner  and  Swasey  Group  composed  of  Harold  Babrov, 
Richard  Tourin,  and  Burton  Krakow  have  made  many  contributions 
to  the  solution  of  the  problem  of  utilizing  spectroscopic 
techniques  for  the  study  of  the  radiant  properties  of  hot 
flowing  gases  (Refs.  60  through  75).  It  is  this  group*s 
development  of  the  modified  Curtis-Godson  approximation 
utilizing  a  random  band  model  with  Lorentz  lines  of  constant 
line  strength  and  line  width  which  provides  the  foundation  for 
the  determination  of  temperature  and  partial  pressure  profiles 
by  the  process  described  in  this  report. 

There  is  also  considerable  interest  in  the  spectral 
radiative  properties  of  carbon  dioxide  as  indicated  by  the 
work  of  Malkmus  in  the  General  Dynamics  Group  (Ref .  49  and  50) . 
Among  others  who  have  been  interested  in  developing  useful 
band  model  techniques  to  molcules  of  the  atmosphere  and 
combustion  are  Oppenheim  and  Ben-Aryeh  (Ref .  76) .  They  have 
been  particularly  concerned  with  low  resolution  measurements 
to  determine  band  model  parameters  for  spectral  emissivity 
calculations  of  the  4.3|j,  band  of  CO2  (Ref.  77  and  78)  .Ellis 
and  Schurin  (Ref .  79)  have  also  made  integrated  intensity 
measurements  on  CO2  bands.  The  near  infrared  is  not  the  only 
spectral  region  for  which  total  band  absorptances  have  been 
measured  and  for  which  band  model  calculations  have  been 
performed.  The  15|i  CO2  band  is  of  considerable  interest  Refs, 
(80  through  84) .  It  was  agreed  that  this  band  would  not  be 
investigated  for  this  project  because  instrumentation  is  not 
available . 

There  are  two  reference  monographs  intended  solely  to 
provide  fundamental  Information  on  the  properties  of  the 
water  vapor  bands  at  2,7\x  (Ref.  85)  and  on  the  1.9ij,  and  6.3iji 
water  vapor  bands  (Ref.  86)  in  terms  of  line  absorptions  and 
integrated  transmittance  at  any  frequency  for  Infinite 
resolution.  Although  the  line  parameters  listed  herein  are 
vital  for  the  interpretation  of  high  resolution  transmittance 
measurements  they  play  no  role  in  the  work  of  the  present 
report  which  is  concerned  only  with  low  resolution  band 
emission-absorption  measurements . 

Tables  of  line  parameters  such  as  the  above  Monographs 
depend  for  their  development  upon  high  resolution  spectroscopic 
data.  Considerable  effort  along  these  lines  has  been  expended 
in  the  past  few  years  by  several  capable  researchers  whose 
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efforts  are  reported  in  the  References  73  and  74  and  85 
through  93.  Although  there  was  Initial  consideration  in  the 
present  project  to  undertake  high  resolution  experiments 
utilizing  an  instrument  such  as  the  Ohio  State  spectrometer 
for  the  purpose  of  preparing  an  atlas  of  heated  water  vapor, 
it  was  later  recognized  that  such  information  was  readily 
available  in  the  open  literature.  Such  high  resolution  spectra 
were  to  be  used  in  selecting  the  most  appropriate  spectral 
region  for  further  investigation  and  for  ascertaining  the 
validity  of  the  various  band  models  for  the  chosen  spectral 
region.  These  things  were  accomplished  on  the  basis  of  other 
considerations  such  as  instrumentation  availability,  difficulty 
and  expense  of  designing  data  acquisition  systems  for  the  test 
cell  environment,  and  the  relative  ease  of  utilizing  low 
resolution  emission-absorption  data  for  band  model  calculations. 
This  appears  to  be  a  reasonable  decision  for  the  present  effort 
but  it  is  not  intended  to  reject  the  usefulness  and  necessity 
of  high  resolution  measurements  for  similar  projects  in  the 
future . 

The  technique  of  total  band  absorptance  and  emittance 
measurements  and  the  utilization  of  band  models  to  predict 
radiance  and  transmittance  are  not  limited  to  H2O  and  CO2  but 
have  been  applied  to  such  molecules  as  HCt  (Ref .  94) ,  CO 
(Refs.  93  and  96),  and  N2O  (Refs.  97  and  98).  The  Elsasser 
strong  line  approximation  was  found  suitable  for  both  CO  and 
N2O. 

Theoretical  examinations  of  the  validity  of  such  concepts 
as  Beer's  Law  (Ref.  99),  the  Doppler  Analog  of  the  Elsasser 
Model,  (Ref.  100),  the  Voight  line  profile  (Ref.  101  and  102), 
the  net  interval  for  calculating  spectra  (Ref.  103),  and  the 
instrumental  effects  of  slit  functions  (Ref.  104),  are  all 
Important  and  continuing  investigations  supporting  the  efforts 
of  experimentalist  to  acquire  meaningful  data  and  hopefully 
helping  everyone  to  a  better  understanding  of  the  physical 
processes  taking  place  in  radlatlng-absorbing  gases. 

The  literature  would  not  be  complete  without  at  least 
one  comparision  between  a  band  model  and  line  by  line 
calculations.  Such  a  comparison  has  been  reported  by  A. 

Goldman  and  T.  G.  Kyle  (Ref.  105)for  the  statistical  band 
model  and  line  by  line  calculation  with  application  to  the 
9.6|i  ozone  and  the  2.7p,  water  vapor  bands.  Their  conclusion 
was  that  the  agreement  between  the  two  types  of  calculation  is 
much  better  than  the  agreement  of  either  of  them  with  experi¬ 
ment  .  Thus  a  refinement  of  the  line  parameters  (high 
resolution  data)  is  needed  in  order  to  obtain  better  agreement 
between  the  calculations  and  the  experimental  results.  In 
general,  however,  the  transmittance  averaged  over  intervals 
greater  than  10  cm-i  for  the  2.7p,  H2O  band  is  represented 
very  well  by  the  statistical  model.  This  supports  and 
correlates  with  the  effort  of  the  present  project. 
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SECTION  III 
BAND  MODEL  THEORY 

The  radiant  energy  emitted  by  a  hot  gas  is  the  integral  of 
contributions  from  all  the  spatial  regions  within  the  volume  of 
the  gas  and  all  the  spectral  intervals  over  which  it  is  active. 
Spatial  integration  by  numerical  methods  is  feasible  because 
spatial  variations  are  usually  gradual  enough  to  allow  utili¬ 
zation  of  large  Intervals.  The  spectral  integration,  however, 
is  often  not  possible.  To  predict  accurately  the  radiation  from 
hot  gases,  the  spectral  absorption  coefficient,  k^,  as  a 
function  of  frequency  must  be  known.  In  molecular  band  radiation 
k^  is  a  rapidly  varying  function  of  frequency,  and  it  is 
essential  for  numerical  integration  that  k^  be  constant  over  each 
frequency  interval.  In  many  situations  the  frequency  interval  is 
so  small  that  it  has  been  desirable  to  develop  techniques  for 
averaging  over  spectral  fine  structure.  Band  models  provide  the 
averaging  processes  by  substituting  for  the  true  spectrum  an 
infinite  array  of  spectral  lines  of  uniform  statistical  properties. 
Four  types  of  band  models  are  particularly  useful  and  will  be 
briefly  discussed:  1)  The  Elsasser  model,  a  band  composed  of 
identical  uniformly  spaced  lines;  2)  The  statistical  model,  a 
band  composed  of  spectral  lines  with  arbitrary  intensity  dis¬ 
tribution  and  with  random  spacing  between  the  lines;  3)  The 
random  Elsasser  model,  a  band  composed  of  several  groups  of 
lines,  each  of  which  forms  an  Elsasser  band,  but  the  groups  are 
superposed  with  random  spacing;  and  4)  The  quasi-random  model, 
a  band  composed  of  spectral  lines  whose  intensities  and  vari¬ 
ation  of  spacing  from  line  to  line  are  accurately  simulated. 

The  band  models  lend  themselves  to  further  approximations 
in  that  there  are  three  limits  in  which  asymptotic  forms  of  the 
equations  defining  each  model  can  be  obtained.  These  are:  1) 

The  weak-line  approximation,  2)  The  strong- line  approximation, 
and  3)  The  non-over lapping-line  approximation. 

Band  models  were  designed  for  and  are  successfully  applied 
to  the  calculation  of  radiation  of  homogeneous  gases.  In  many 
cases,  however,  the  hot  gases  are  no  longer  homogeneous  but 
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have  gradients  of  temperature,  pressure,  and  concentration.  The 
Curtls-Godson  approximation  provides  a  technique  for  calculating 
the  radiance  and  transmittance  of  hot  Inhomogeneous  gases  by 
determining  the  necessary  spectral  parameters  solely  from 
homogeneous  gas  data.  In  the  Curtls-Godson  approximation  the 
transmittance  of  the  Inhomogeneous  path  is  obtained  by  substi¬ 
tuting  a  hypothetical  homogeneous  path  with  the  same  trans¬ 
mittance.  The  homogeneous  path  can  then  be  treated  by  establish¬ 
ed  band  model  methods. 

The  Single  Spectral  Line 

This  may  be  thought  of  as  the  most  simple .band  model.  It 
is  the  situation  in  which  the  spectral  lines  are  so  narrow  and 
far  apart  that  there  is  effectively  no  overlapping  of  lines  and 
each  line  can  be  treated  as  an  entity.  This  situation  will  be 
discussed  In  detail  In  order  to  introduce  concepts,  define  para¬ 
meters,  and  discuss  limitations  that  are  fundamental  to  band 
model  theory. 

In  view  of  the  complicated  conditions  that  are  often  pre¬ 
sent  In  a  source  of  radiation  It  is  more  fruitful  to  discuss 
spectral  lines  In  terms  of  absorption  rather  than  emission  since 
conditions  inside  an  absorption  cell  can  be  kept  relatively 
simple  and  controllable.  If  parallel  light  from  a  source 
emitting  a  continuous  spectrum  of  intensity  I^.^^  is  sent  through 
an  absorption  cell  containing  a  monoatomlc  gas,  the  intensity  of 
the  transmitted  light,  I.  .  may  show  a  frequency  distribution 
similar  to  that  of  Figure  1.  When  this  is  the  case  the  gas  is 
said  to  possess  an  absorption  line  at  the  frequency  Note 

that  the  spectral  line  has  a  definite  width  or  spread  of 
frequencies  around 

The  absorption  coefficient,  k.^,  of  the  gas  is  defined  by 
the  Lambert -Beer  relation 

dl  =-k^Idt 

which  is  the  flux  absorbed  frcnn  a  beam  of  radiation  of  frequency 
u  of  unit  cross  section  as  it  traverses  a  path  length  d^  through 
the  gas.  When  1.  is  measured  in  centimeters,  k.^^  is  expressed 
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in  reciprocal  centimeters.  The  spectral  transmittance,  is 

defined  as  the  ratio  of  the  transmitted  flux  to  the  incident 
flux  at  frequency  u  and  may  be  written  in  terms  of  the  absorp¬ 


tion  coefficient  as 


h/i\  = 


-k  -t 

e  x> 


From  this  equation  and  Figure  1  we  may  obtain  as  a  function 
of  frequency 

^  -  T  h/i°  (2) 


and  plot  it  as  in  Figure  2.  The  halfwidth  of  the  line,  7,  is 
shown  in  Figure  2  and  is  defined  as  the  total  width  of  the 
absorption  coefficient  curve  at  the  point  where  the  absorption 
coefficient  is  one-half  its  maximum  value.  In  general  the 
absorption  coefficient  of  a  gas  is  given  by  an  expression  involv¬ 
ing  a  function  of  the  frequency  u,  a  definite  value  of  k„„„, 

xZiElX 

and  the  halfwldth  7  all  of  which  depend  on  the  nature  of  the 
molecules  of  the  absorbing  gas,  their  motion,  and  their  inter¬ 
action  with  one  another  or  with  foreign  molecules. 

The  intensity  of  a  single  line  can  be  measured  in  absorption 
by  determining  the  Integrated  absorption,  S,  defined  as  the 

integral  of  the  absorption  coefficient, 

00 

S  =  /  (3) 


which  is  simply  the  area  under  the  absorption  coefficient  curve, 
usually  given  in  units  of  cm”^.  jn-  general  Beer’s  law  does 
not  hold,  i.e.,  the  absorption  is  not  proportional  to  concen¬ 
tration.  However  for  gas  at  densities  sufficiently  low  that 
the  mechanism  for  absorption  is  not  altered,  the  integrated 
absorption  is  directly  proportional  to  the  n\imber  of  molecules 
of  the  absorber  in  the  optical  path;  thus, 

S  =»  S°P  (4) 

where  S°  is  the  line  strength  in  units  of  cra~®atm~^  and  p  is 
the  pressure  in  atm  of  the  absorber.  The  line  strength  is 
related  directly  to  the  theoretical  transition  probability. 
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The  spectral  absorptance,  a^,  which  provides  a  useful  way 
of  expressing  the  absorption  of  radiation  at  frequency  o,  is 
defined  as  the  ratio  of  the  difference  between  the  incident 
'flux  and  the  transmitted  flux  to  the  incident  flux  at  fre¬ 
quency  o .  Thus 


a 


V 


(5) 


•0 

When  conditions  are  such  that  local  ‘thermodynamic  equilibrium 
can  be  assumed,  the  spectral  absorptance,  a^,  may  be  equated 
to  the  spectral  emissivity,  by  Kirchhoff's  Law.  It  is 

this  assumption  which  makes  possible  the  calculation  of  the 
radiant  energy  emitted  by  a  hot  gas  utilizing  the  concepts  of 
absorption  spectroscopy. 

The  integrated  absorptance.  A,  over  a  finite  wavenumber 
interval,  Au,  is  defined  as 


A  '  (6) 


The  integrated  absorptance  is  also  called  the  equivalent  width, 

W.  is  contributed  to  by  one  line  only,  then  the  integral 

of  the  spectral  absorptance  over  all  frequencies  is  finite  and 
the  equivalent  width  is  given  by 

W  =  7  (1  -  (7) 

o 

which  is  simply  the  integral  of  the  fractional  flux  removed  from 
the  Incident  flux,  i.e.,  the  area  under  the  absorption  curve. 

The  equivalent  width  is  a  function  of  the  amount  of  absorber 
and  the  relationship  between  the  two  is  called  the  curve  of 
growth. 

Either  the  spectral  absorptance  or  the  spectral  transmittance 
is  the  quantity  ideally  sought  in  single  line  absorption 
experiments.  From  the  resulting  absorption  coefficient  curve 
the  Integrated  absorption  S  may  then  be  measured  and  compared 
directly  with  theory.  These  preferred  measurements, which 
require  the  techniques  of  high  resolution  spectroscopy,  are 
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very  difficult  partly  because  of  undetermined  effects  of  the 
instrument  on  the  spectral  line  shapes.  Consequently  it  is 
usually  the  equivalent  width,  W,  which  is  actually  measured 
because  experimental  conditions  can  in  most  cases  be  so  chosen 
that  the  equivalent  width  is  independent  of  the  slit  function 
of  the  spectrometer.  The  equivalent  width  is  then  related  to 
the  integrated  absorption  (and  hence  the  theory)  through  an 
absorption  coefficient  curve.  These  curves  have  a  shape  and 
a  halfwidth  for  each  individual  line  which  is  determined  by 
the  interactions  of  the  molecules  in  the  gas.  Thus  the  line 
shape  and  halfwidth  are  dependent  on  the  temperature  and 
pressure,  i.e.,  upon  the  mechanisms  responsible  for  broadening 
the  line.  These  line  broadening  mechanisms  can  be  categorized 
as  natural  broadening,  Doppler  broadening,  and  collision 
broadening.  The  natural  width  is  inversely  proportional  to 
the  finite  lifetime  of  the  excited  states  and  compared  to  other 
processes  can  be  neglected  for  the  measurements  in  the  infrared 
spectral  region.  For  low  pressure  situations  the  principal 
line  shape  factor  is  the  Doppler  effect  where  the  absorption 
coefficient  is  given  by 

k  =  ( g  ^  (8) 

u,D  "  it"  7d 

and  the  Doppler  halfwidth  is  given  by 


2u  j.2WI-ln2-,h  /qn 

'>'d  =  tr  [— H — ’ 

where  M  is  the  molecular  weight 

As  the  pressure  is  increased  the  principal  process  respon¬ 
sible  for  determining  the  line  shape  is  that  of  binary  colli- 

» 

sions.  In  this  pressure  regime  a  satisfactory  representation 
for  the  line  shape  is  given  by  the  Lorentz  line  shape  formula 
which  in  the  infrared  is 


(10) 
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where  7^  is  the  Lorentz  halfwldth.  The  Lorentz  halfwldth  is 
proportional  to  the  nvunber  of  collisions  per  second  interrupting 
the  radiation  absorption  process  and  so  formulas  can  be  obtained 
from  kinetic  theory  for  this  halfwldth.  To  a  good  approximation 
one  can  simply  say  that  7^  is  proportional  to  the  total  pressure 
of  a  pure  gas,  and  write 

7°?  (11) 

Expressions  are  also  available  for  lines  which  are  broadened  by 
a  combination  of  Doppler  and  collision  broadening. 

These  are  variations  of  the  Voigt  Formula  which  is  given  by 


-y 


-00  a  +(w  -  y)‘ 


dy 


(12) 


where 


y  =  ^  y/US  ~  an  integration  interval 

k  =  maximum  absorption  coefficient 
when  only  Doppler  Broadening  is 
present 

a  -N,  —  v/^n2 
"■  '^D 

The  mixed  Doppler  and  Lorentz  line  shape  is  not  amenable  to  an 
analytical  solution  but  is  solved  by  numerical  techniques. 

Three  different  theoretical  expressions  for  the  absorption 
coefficient,  k^,  have  been  given,  each  representing  fairly 
accurately  the  shape  of  the  spectral  line  if  the  appropriate 
conditions  are  satisfied.  If  the  shape  and  strength  of  a  line 
are  known,  its  equivalent  width  may  be  calculated  from  Equation 
7.  For  lines  with  the  Lorentz  line  shape  the  equivalent  width 
is  given  by  the  Ladenburg  -  Heiche  equation 

W  =  2ir7Lf(x)  (13.) 
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The  function  f (x)  is  defined  by  an  expression^  in  Bessel  functions 
as 

f(x)  =  xe“^[J^(ix)  -  i  Jj(ix)]  (14) 

where 

X  =  S^/2ir7j^  =  S°V2t:7°  (15) 

For  a  pure  gas  x  is  independent  of  pressure  and  can  be  varied 
only  by  varying  the  path  length,  t.  For  a  gas  mixture  x  can 
be  varied  for  a  particular  line  at  a  constant  path  length  by 
varying  the  partial  pressure  of  the  absorbing  gas.  The  import¬ 
ance  of  the  variation  of  x  becomes  evident  when  it  is  recognized 
that  in  order  to  determine  both  the  line  strength  S°,  and  the 
halfwldth  7°,  for  a  particular  spectral  line  from  the  measure¬ 
ment  of  equivalent  width,  it  is  necessary  to  determine  W  for 
more  than  one  value  of  x.  This  becomes  apparent  when  the  limit¬ 
ing  expressions  for  the  equivalent  width  are  examined.  For 
small  values  of  x 

f (x)  ^  X  (1  -  x/2  +  .  .  .)  (16) 

and  the  linear  region  of  the  curve  of  growth  is  approached  where 
the  equivalent  width  is  given  by 

Vf®  =  I  ^  S°t[l  -  W®/4Tr7®  +  .  .  .]  (17) 

For  large  values  of  x  the  square  root  region  is  approached  where 

f(x)  ^  (f^)2  [1  +  l/8x  +  .  .  .]  (18) 

and 


W°  ^  2(S°  Tt)2[l  -  ir(7%°)2+  .  .  .  ]  (19) 

Thus  there  are  two  asymptotic  approximations  for  the  equivalent 
widths  of  isolated  spectral  lines: 

Linear  Approximation 

=  S°l  X  «  1  (20) 
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Square  Root  Approximation 

=  2(S°7°^)2  X  >  3  (21) 

The  determination  of  the  spectral  parameters  of  line  strength 
and  line  width  is  important  not  only  to  the  theory  of  isolated 
spectral  lines  but  also  to  band  model  theory  because  the  band 
model  parameters  may  be  determined  from  these  more  fundamental 
molecular  quantities. 

Approximations  to  Band  Models 

Before  discussing  the  individual  band  models,  a  brief 
preview  of  the  approximations  to  the  models  will  be  given 
because  each  model  will  be  considered  in  the  light  of  these 
general  limiting  cases. 

Dr.  G.  N.  Plass  proposed  three  approximations  (Ref.  JOSA 
50,  869  (I960)) that  are  Intended  to  permit  predictions  of  band 
absorptance  for  virtually  any  range  of  pressure  and  concentration 
of  absorbing  gas  in  terms  of  relatively  simple  mathematical 
expressions  involving  the  parameters  X  and  p.  X  =  Sii/2ir7,  is  a 
measure  of  the  ratio  of  the  line  strength,  S,  to  the  line  half- 
width,  7,  for  an  absorber  concentration,  ix.  p  *=  2ir7/d,  is  the 
ratio  of  the  line  halfwidth,  7,  to  the  mean  line  spacing,  d. 

The  simple  expressions  for  the  band  models  apply  to  limiting 
curves  on  logarithmic  plots  giving  absorptance  as  a  function  of 
X,  or  p,  or  various  products  of  these  two  parameters.  The  three 
Plass  approximations  are  the  strong  line,  the  weak  line,  and  the 
nonoverlapping  line  approximation.  Each  of  these  approximations 
has  a  well  defined  region  of  validity  discussed  briefly  below 
and  curves  are  presented  in  the  original  paper  which  show  these 
regions  of  validity  for  the  different  band  models. 

1 .  Strong  Line  Approximation 

When  the  absorptance  is  virtually  complete  near  the  centers 

of  the  strongest  lines  in  the  band,  the  band  absorptance  is  a 

2 

function  of  the  single  variable  p  X  where 

p^x  =  (^)^  (|^^)  =  2177  S^/d^  (22) 
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The  strong  line  approximation  Is  defined  to  be  this  limiting 
absorption  curve.  This  approximation  is  valid  even  when  the 
spectral  lines  overlap  and  should  not  be  confused  with  the  square 
root  approximation  which  is  valid  only  when  the  spectral  lines 
do  not  overlap.  The  absorption  as  calculated  from  the  strong 
line  approximation  depends  on  the  arrangement  of  the  spectral 
lines  in  the  band.  It  is  different  for  a  random  spacing  of  the 
lines  than  for  a  regular  spacing  and  so  the  shape  of  the 
experimental  curve  should  serve  as  a  guide  to  the  selection  of 
an  appropriate  band  model.  The  strong  line  approximation  is 
useful  in  extrapolating  laboratory  data  to  large  values  of 
absorber  concentration ,  \i,  and  small  values  of  pressure,  p. 
Laboratory  measurements  in  this  approximation  would  be  made  in  a 
long  absorption  cell  at  low  pressure. 

2 .  The  Weak  Line  Approximation 

When  the  spectral  absorptance  is  small  at  all  frequencies 
in  a  band  including  those  near  the  centers  of  the  strongest 
lines,  the  absorptance  is  a  function  of  the  single  variable  px 
where , 

PX  =  (  )  =  SM,/d  (23) 

This  limiting  absorption  curve  is  called  the  weak  line  approxi¬ 
mation  and  is  valid  even  when  the  spectral  lines  overlap.  Thus 
it  should  not  be  confused  with  the  linear  approximation  which 
is  valid  only  when  there  is  no  overlapping  of  the  spectral  lines. 
This  approximation  is  Independent  of  the  particular  arrangement 
of  the  spectral  lines  in  a  band  and  so  is  nearly  the  same  for 
all  models.  The  weak  line  approximation  is  useful  in  extrapol¬ 
ating  to  small  values  of  absorber  concentration  and  large 
pressures  and  laboratory  measurements  would  be  made  in  a  very 
short  cell. 

3 .  The  Nonoverlapping  Line  Approximation 

When  the  value  of  3,  which  is  proportional  to  the  pressure, 
is  sufficiently  small  that  there  is  no  overlapping  of  lines  the 
band  absorption  is  simply  the  sum  of  the  absorptances  of 
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Individual  lines  as  given  by  the  Ladenburg-Reiche  expression. 

The  only  requirement  for  the  validity  of  this  approximation  is 
that  the  spectral  lines  do  not  overlap  appreciably.  The  non¬ 
overlapping  approximation  is  the  limiting  curve  in  a  logarithmic 
plot  of  A/P  for  constant  p  values  vs.  X.  This  curve  has  a 
linear  region  for  weak  nonoverlapping  lines  and  a  square  root 
region  for  strong  nonoverlapping  lines.  The  nonoverlapping 
approximation  is  useful  in  extrapolations  to  small  p.  and  small 
P. 

Plass  has  shown  that  the  approximations  should  be  valid 
over  practically  all  ranges  of  X  and  p  for  the  Elsasser  model 
and  over  most  ranges  of  X  and  p  for  the  statistical  model  and 
that  predictions  of  absorption  to  within  ±  10%  are  possible. 
Plass  also  describes  laboratory  studies  appropriate  to  the 
selection  of  the  correct  approximation  for  the  particular  situ¬ 
ation  of  interest. 

The  Elsasser  Model 

The  Elsasser  band  model  consists  of  an  infinite  array  of 
spectral  lines  of  equal  intensity  and  equal  halfwidths  spaced 
at  equal  intervals.  It  has  been  shown  to  be  most  suitable  for 
linear  molecules.  This  model  is  based  on  the  Lorentz  line  shape 
and  is  indeed  simply  an  infinite  repetition  of  a  Lorentz  line 
so  that  the  absorption  coefficient  is  given  by 


•0 


00 

=  z 

n=-oo 


(n  -  nd)'^  +  7_ 


(24) 


where  d  is  the  line  spacing 

It  is  possible  to  express  this  infinite  sum  in  terms  of 
periodic  and  hyperbolic  functions, 


where 


k 


_  s  sinh  p 
u  H  cosh  S  -  COS  z 


3  =  2ir7j^/d 


z  =  2ir'o/d 


(25) 
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If  the  averaging  interval  A'o  is  taken  as  one  period  of  the  band 
(An  =  2ir)  the  exact  expression  for  the  integrated  absorptance 
or  equivalent  width  of  an  Elsasser  band  is 

A  =  ^  /  (1  -  e"^n'^)dz 

—IT 


IT 

=  1  -  i-  /  exp[  (-6x)  - sinhg -  i^z  (26) 

2¥  .  i-L  V  /  cos 6  -  cos  Z  ■■ 

“IT 

where  px  =  st/d 

This  integral  cannot  be  evaluated  in  closed  form  and  must  be 
integrated  numerically.  Several  solutions  based  on  expansions 
of  the  integrand  have  been  developed;  however,  it  is  often  quite 
useful  to  work  with  the  asymptotic  approximations  which  are 
valid  for  certain  conditions. 

The  weak  line  approximation  is  the  asymptotic  limit  as 
P  00.  Then  sinhp  ->  oo  and  the  integrated  absorptance  becomes 

A  =  1  -  e-P*  -  1  -  (27) 

Since  the  parameter  p  measures  the  ratio  of  the  line  width  to 
the  distance  of  the  neighboring  line',  in  this  limit  lines  over¬ 
lap  and  there  is  no  fine  structure.  The  transmission  is  now 
independent  of  pressure  (Beer^s  Law).  This  approximation  is 
good  whenever  the  absorption  at  the  line  centers  is  small,  i.e., 
when  the  ratio  of  the  line  strength  to  the  line  halfwidth  is 
small  and  for  this  reason  is  referred  to  as  the  weak  line 
appr oximat ion . 

The  strong  line  approximation  is  the  asymptotic  limit  to 
equation  26  when  x>  1.  Under  this  condition  the  absorption  at 
the  line  center  is  usually  complete,  the  halfwidths  are  narrow 
and  the  lines  do  not  overlap  strongly  (p  small) .  The  integral 
becomes , 

A  =  erf(|p^x)2  (28) 

where  the  error  function  is  given  as, 
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9  *' 

erf(t)  =  —  /  e  ^  dt 
aA"  o 

This  approximation  is  particularly  useful  for  long  path  lengths 
and  small  values  of  pressure.  It  differs  from  the  square  root 
approximation  in  that  here  the  lines  may  overlap. 

The  Statistical  Model 

The  statistical  or  Mayer-Goody  model  of  a  band  assumes  a 
random  spacing  of  spectral  lines  and  an  Intensity  that  varies 
from  line  to  line  in  any  manner  which  can  be  represented  by  a 
distribution  function.  An  expression  for  the  integrated  absorp- 
tance  for  this  model  can  be  developed  from  the  following  con¬ 
siderations  . 

Let  A13  be  a  spectral  interval  in  which  there  are  n  lines 
of  mean  distance  d, 

Ai3  =  nd 

where  it  is  assumed  that  any  line  has  equal  probability  of  being 
anywhere  in  the  interval. 

Let  P(S^)  be  the  probability  that  the  i  line  will  have 
an  intensity  between  and  S^+dS^.  The  probability  is  normal¬ 
ized  so  that  , 

00 

f  P(S^)dS^  =  1 
o 

The  mean  transmittance  is  found  by  ayeraging  over  all  positions 
and  all  intensities  of  the  lines,  thus 


T  *  — •  •  •  / 

^  AB  “j 

But  since  all  the  integrals  are  alike  we  may  write 

T  *  P(s)e"^v^ds]“ 


which  may  be  rewritten  as 
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-  -^/^v/PCsXl  -  e-V)ds]" 
or  interchanging  the  order  of  integration 

T  =  [1  -  /P(s)ds  /  (1  -  e~\^)dx>] 

The  integral  over  the  frequency  in  this  equation  is  equal  to 
the  equivalent  width  of  a  single  line  evaluated  over  the  fre¬ 
quency  Interval  Au.  Thus  the  transmittance  can  be  written 

The  average  value  of  the  equivalent  width  over  the  distribution 
of  line  strengths  is 

7)P(S)dS 

j  »  O’ 

where  indicates  some  mean  line  intensity  that  occurs  in  the 
distribution  function  P(S) .  Since  Au  =  nd  the  transmittance 
becomes 

<29) 

If  the  number  of  lines  in  the  band  is  allowed  to  approach  infin¬ 
ity  while  the  mean  spacing  is  held  constant,  this  last  express¬ 
ion  approaches  an  exponential 

T  =  exp[-  p/d]  (30) 

In  terras  of  average  integrated  absorptance  for  the  band  equations 
(29)  and  (30)  become  respectively 

the  absorption  from  a  finite  number  of  spectral  lines,  n;  and 

A  =  1  -  exp[  -  jj/d]  (32) 

the  absorption  from  an  infinite  number  of  randomly  spaced 
spectral  lines. 

The  equations  are  often  examined  for  two  special  intensity 
distributions;  lines  of  equal  intensity  and  line  strength  by 
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Poisson  distribution. 

Where  all  the  lines  have  equal  intensity  equation  (21) 
reduces  to 

A  =  1  -  exp  I  -  pxe“*[J^(ix)  -  i  J^Cix)]  [  (33) 

If  each  of  the  lines  absorbs  weakly  so  that  x  is  small  then  we 
have  the  weak  line  approximation  to  the  statistical  model  with 
all  lines  equally  intense  given  by 

A  =  1  -  e"^^  (34) 

If  the  lines  absorb  strongly  then  we  have  the  strong  line 
approximation  to  the  statistical  model  with  all  lines  equally 
intense  given  by 

A  =  1  -  exp[-(^  P^x)2] 

Many  spectral  bands  may  be  simulated  by  assuming  an 
exponential  distribution  of  line  intensities,  namely, 


P(S)  =  S^’^expC-  S/Sq)  (35) 

With  this  distribution,  Goody  developed  the  following  formula 
for  the  Integrated  absorptance 

A  =  1  -  exp[  -  Px^Cl  +  (36) 

where  =  SQ-t/27r7 

The  weak  line  approximation  to  equation  (36)  is  obtained  when 

Xo  «  1 


A  =  1  -  exp(-  8X^) 


(37) 


The  strong  line  approximation  to  equation  (36)  is  obtained  when 


Xo»l 


A  =  1  -  exp[-(^P^X)2] 


(38) 


\ 
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The  Random  Elsasser  Model 

The  predicted  absorption  for  the  Elsasser  and  statistical 
models  agree  at  small  values  of  x  because  the  absorption  is 
determined  by  the  total  strength  of  all  the  spectral  lines. 
However,  as  x  Increases  the  results  calculated  from  these  two 
models  begin  to  diverge,  with  the  Elsasser  theory  always 
predicting  more  absorption  than  the  statistical  model  for  a 
given  value  of  x.  Thus  another  model  which  is  a  generalization 
of  these  two  models  is  suggested.  The  spectral  lines  of  an 
actual  absorption  band  are  not  arranged  either  completely  at 
random  or  at  regular  intervals.  The  actual  line  pattern  is 
formed  by  the  superposition  of  many  systems  of  lines. 

The  random  Elsasser  band  model  simulates  this  actual 
pattern  by  assuming  that  the  absorption  may  be  characterized 
by  the  random  superposition  of  Elsasser  bands,  each  of  which 
may  have  a  different  line  intensity,  halfwidth,  and  spacing 
between  the  lines.  As  many  different  Elsasser  bands  as  desired 
may  be  combined  by  multiplication  in  this  model  and  thus  the 
weak  spectral  lines  that  contribute  to  the  absorption  can  be 
included  in  the  calculations.' 

The  absorptance  of  N  randomly  superimposed  Elsasser  bands 
is 

N 

A  =  1  -  n  [1  -  Aj.  i(x^,p.)]  (39) 

where  A„  .  is  the  absorptance  of  an  Elsasser  band  with  a  half- 
* ,  1 

width  of  7^,  a  line  spacing  d^,  and  a  line  strength  so  that 

Xi  =  Sj^-t/2ir7^  and  =«  27r7j^/d^ 

The  weak- line  approximation  for  the  absorptance  in  the  random 
Elsasser  model  is  given  by 

N 

A  =  1  -  IT  exp[  -  P.x.  ]  (40) 

i=l  ^  ^ 

The  strong-line  approximation  is 

A  =  1  -  .n^[l  -  erf  (Ip.^x.)^]  (41) 
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As  the  number  of  randomly  superimposed  Elsasser  bands  increases 
the  absorptance  approaches  that  of  the  statistical  model. 

The  Quasi-Random  Model 

The  quasi-random  model  does  not  require  that  the  spectral 
lines  be  either  uniformly  nor  randomly  spaced  but  is 
sufficiently  flexible  to  simulate  the  complex  spacing  and 
intensity  distributions  of  lines  which  are  actually  observed 
in  practice.  It  provides  for  calculating  the  effect  of  wing 
absorption  from  spectral  lines  whose  centers  are  outside  the 
frequency  interval  of  interest,  as  well  as  including  the  con¬ 
tributions  of  the  weak  lines  in  the  interval.  This  is  the  most 
accurate  of  the  band  models  but  it  may  approach  line  by  line 
calculations  in  complexity.  In  this  model  the  general  trans¬ 
mittance  cannot  be  expressed  in  terms  of  elementary  functions 
and  an  a  priori  knowledge  of  the  band  structure  is  required. 
However,  because  of  its  greater  sophistication  this  band  model 
provides  the  most  realistic  but  perhaps  unnecessarily  complicated 
treatment  for  calculating  the  temperature  and  concentration 
profiles  of  hot  gases  from  measurements  of  the  radiant  energy. 

In  the  quasi-random  model  the  absorptance  is  calculated 
first  for  a  frequency  interval  that  is  much  smaller  than  the 
Interval  of  interest.  This  localizes  the  stronger  lines  to  a 
narrow  interval  around  their  actual  positions  and  prevents  the 
introduction  of  spurious  overlapping  effects.  The  absorptance 
for  each  of  the  N  spectral  lines  in  the  interval  is  calculated 
separately  and  the  results  combined  by  assuming  a  random 
placing  of  the  lines  within  the  small  Interval.  The  absorptance 
values  for  all  of  the  small  intervals  that  fill  the  large 
interval  of  Interest  are  averaged  to  obtain  the  final  value  for 
the  absorptance.  The  absorptance  for  the  quasi-random  model  is 

A  =  i  2  A.  (42) 

^  j=l  ^ 

where  Aj  is  the  absorptance  of  each  of  the  L  smaller  intervals 
and  is  calculated  from 
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Where  ^  is  the  single  line  absorptance  over  the  finite 
interval  D  and  N  is  the  number  of  lines  in  the  frequency  inter¬ 
val  j  . 

The  weak-line  approximation  for  the  quasi-random  model  is 
given  by 

A  =  i  2  |l  -  n  Cl  -  e“  Pi^i)[  (44) 

^  j=l  «  i=l  ’ 

the  strong  line  approximation  is, 

A=l2  L-  -  erf(Ci)l 

^  J=1  ' 

Where 

C 

The  quasi-random  model  is  developed  and  discussed  in  detail  by 
Wyatt,  Stull,  and  Plass  (Ref.  17). 

The  Curtis  Godson  Approximation 

Band  models  have  been  applied  with  considerable  success 
to  calculate  the  transmittance  of  homogeneous  gas  sample.  In 
many  practical  situations,  however,  hot  gas  specimens  have 
gradients  of  temperature,  pressure  and  concentration.  It  then 
becomes  necessary  to  extend  the  band  model  considerations  to 
the  calculation  of  the  transmittances  of  Inhomogeneous  gases . 
There  are  a  number  of  possible  ways  of  developing  the  required 
approximation  but  perhaps  the  most  useful  is  the  Curtls-Godson 
approximation.  (Ref .4)  In  the  Curtis-Godson  approximation, 
the  non-homogeneous  path  is  approximated  by  an  homogeous  path 
by  assigning  an  effective  amount  of  absorbing  gas  and  an 
effective  pressure  referred  to  a  standard  temperature  so  that 
the  two  paths  will  absorb  exactly  the  same.  This  implies  that 
a  single  line  absorbing  along  a  non-homogeneous  path  can  be 
simulated  by  a  line,  whose  intensity  and  width  can  be  adjusted. 
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absorbing  along  a  homogeneous  path.  The  adjustment  of  intensity 
and  width  is  made  to  give  the  correct  absorption  in  the  strong 
and  the  weak  line  regions.  Since  the  interest  here  is  in  the 
application  to  combustion  phenomena  the  development  followed 
will  be  that  of  the  Warner  and  Swasey  group  who  utilized  a 
finite  difference  scheme  to  modify  the  Curtis-Godson  method. 

This  simplification  allowed  the  use  of  inhomogeneous  paths 
that  were  composed  of  a  series  of  smaller  homogeneous  zones. 

The  optical  path  is  divided  into  r  isothermal  zones  as 
shown  in  Figure  3.  The  spectral  radiance  of  this  inhomogeneous 
column  of  gas  is  then  given  by 


n 


2  -1  —1 

in  units  of  W[cm  cm  str] 


(46) 


where 

i  =  labels  an  Isothermal  zone  in  the  line  of 

sight;  i  increases  with  the  distance  of  the 
zone  from  the  detector. 

J  =  labels  a  spectral  Interval. 

Uj  =  frequency  of  the  midpoint  of  the  spectral  interval 

Au .  =  frequency  spread  of  the  spectral  interval 
3 

Ny^(u.  T.  )  =  spectral  radiance  of  a  Black  Body  at  the  frequency 
D  j »  1 

n,  and  temperature  T.  (Planck's  Function) 

J  ^ 

€ , (u .  Au.)  =  spectral  emissivity  at  the  freq.  u.  of  zone  i. 
t  3 »  J  3 

Under  conditions  such  that  local  thermodynamic  equilibrium 

can  be  assumed,  the  spectral  emittance  may  be  equated  to  the 

spectral  absorptance  by  Kirchhoff's  Law 


6.  (t> .  AC  )  =  a.  (u .  Ac.)  -  a,  ^  (u  Au  ) 
iJ^J  ^  J  f  J 


(47) 


or  since  =  1  - 


6.  (t> .  At>.)  =  T.  ,  (u  .  AO.)  -  T.  (u.,A^3.) 
J  >  J  ±"-i  J  f  J  J  J 


(48) 
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Here  t.(u.  Au.)  is  the  transmittance  of  the  column  of  gas  from 
^  3  f  3 

the  beginning  of  the  first  zone  (the  zone  nearest  the  detector) 
to  the  end  of  zone  i.  Similarly  j^C^dfAUj)  is  the  trans¬ 
mittance  from  the  beginning  of  the  first  zone  to  the  end  of 
zone  i  -  1.  taken  as  unity. 

Rewriting  equation  46,  the  spectral  radiance  of  the  inhomogeneous 
column  of  gas  is 


n 

2 

i=l 


T 


T.  (V).  Ao,)]  (49) 

i  J  I  J 


The  radiant  flux,  or  the  irradiance,  is  the  power  per  unit  area 
at  an  element  of  the  detector  and  is  given  by 


m 

H  =  /  2  N(u.  Au  .  )A'o  .cos0dw  (50) 

w  3=1  J  J 

-2 

in  units  of  Wcm 

These  'radiance  equations  provide  a  satisfactory  means  of 

developing  the  calculation  whenever  it  is  possible  to  identify 

isothermal  zones.  The  integral  over  the  solid  angle  can  be 

carried  out  by  numerical  integration.  The  zonal  transmittances , 

T.(t),  At).),  are  the  unknowns  and  are  to  be  determined  by 
1  J »  3 

utilizing  a  suitable  band  model  together  with  the  Curtis-Godson 

approximation.  Each  line  in  the  spectral  interval  Au  -  is 

3 

characterized  by  its  position  v  its  strength  S  and  its  halfwidth 

y.  The  band  model  parameters  are  the  effective  line  strength 

S/d  and  the  effective  halfwldth  y/d  where  d  is  the  frequency 

spacing  between  lines  for  the  interval  Au . . 

3 

In  applying  the  Curtis-Godson  approximation  one  first 
determines  the  values  of  the  band  model  parameters  (S/d)^ 

and  (7/d)j^  for  each  of  the  isothermal  zones  constituting  the 
Inhomogeneous  sample.  The  subscript  h  identifies  any  individual 
zone.  The  band  model  parameters  of  the  equivalent  hypothetical 
homogeneous  path  are  then  defined  in  terms  of  the  parameters 
of  the  individual  zones  by  the  equations 


27 


AEDC-TR-71-78 


(S/d)l  =  2  (S/d).  (51) 

h  "  " 

(S/d) t (7/d)  =  2  (S/d)jj'tjj(7/d)ij  (52) 

where  t  is  the  length  of  an  optical  path. 

It  is  important  to  observe  that  there  is  a  different 
inhomogeneous  path  through  the  gas  sample  for  each  transmittance 
T.(u.  Au.).  This  is  because  t.(u.-Au.)  is  defined  as  the 
transmittance  of  the  inhomogeneous  column  of  gas  from  the 
beginning  of  the  first  zone  to  the  end  of  zone  i.  Hence  there 
must  be  an  equivalent  homogeneous  path  hypothesized  for  each 
transmittance  t . (u  .  Au.)  and  its  band  model  parameters  are 

1  j  j  j 

defined  by  the  equations, 

i 

[(S/d)^].  =  2  (S/d).^.  (53) 

1  h=l  “  “ 

i 

[(S/d)>t(7/d)]^  =  2  (S/d)  (7/d)  (54) 

h— 1 

The  transmittance  through  a  given  non-homogeneous  sample 
is  calculated  by  choosing  a  particular  band  model.  This  may  be 
done  on  the  basis  of  the  known  character  of  the  spectrum  of  the 
gas  under  examination.  For  gases  of  particular  interest  in 
this  study,  namely,  hot  water  vapor  and  carbon  dioxide  it  has 
been  shown  that  transmlttances  can  be  predicted  satisfactorily 
using  the  statistical  band  model  with  constant  line  widths . 
According  to  this  band  model,  the  transmittance  of  a  homogeneous 
gas  over  a  spectral  interval  containing  many  lines  can  be 
expressed  as 

-^n  T  =  2-ir(7/d)f  (x)  (55) 

where  x  is  given  by 

X  =  (S/d)-t/2ir(7/d)  (56) 

The  form  of  f(x)  depends  on  the  probability  distribution  function 
used  for  the  line  strengths .  f (x)  is  the  Ladenburg-Reiche 
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function  for  a  delta  function  distribution  where  all  the  lines 

are  assumed  to  have  equal  intensity.  The  present  development 

will  be  based  on  this  type  of  line  strength  distribution. 

To  determine  the  transmittance  t.(u.  A^3.)  of  an  inhomogen- 

^  J  >  j 

ecus  column  of  gas  consisting  of  i  isothermal  zones  one  calculates 
the  transmittance  of  the  equivalent  hypothetical  homogeneous 
column  whose  band  model  parameters  are  given  by  equations  53 
and  54 .  Thus , 


-^n  AUj)  =  2Tr(7/d)^  f(x^) 


(57) 


The  expression  for  (7/d)^  in  terms  of  the  band  model 
parameters  of  the  isothermal  zones  is  obtained  by  diving 
equation  54  by  equation  53 


and  x^  is 


[  (S/d) ^ (7/d)] 
^  [(S/d)^]. 

[(S/d)^]^ 

^  [  2iT7/d]  ^ 


(58) 


(59) 


Substituting,  we  have  the  transmittance  in  terms  of  zonal  band 
model  parameters. 


-tn  T.  (u.  Au  .)  =  2ir 

^  J  I  J 


Z  (s/d)^t^(7/d)^ 
h  =1 


2  (S/d)  I 

h=l  “  “ 


ri=*l 


(60) 


We  now  return  to  Ladenburg-Reiche  expressions  to  obtain  the  band 
model  parameters  of  the  isothermal  zones  in  terms  of  trans- 
roittances  within  those  individual  zones.  Thus, 


(7/d)h 


--tn 

i^irf  (Xjj) 


(61) 
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and 


(62) 


The  Transmittance  t-Cij.,  A^3.)  of  the  Inhomogeneous  path  can  now 

1  j  j 

be  determined  in  terms  of  the  zonal  transmittances  and  zonal  x's 


This  equation  can  be  used  for  any  value  of  x  but  x  must  be  known, 
f (x)  is  tabulated  for  numerical  values  of  its  argument. 

Since  there  are  two  asymptotic  limits  on  the  Ladenburg- 
Reiche  function  the  above  equation  reduces  to  forms  in  which 
is  independent  of  Xj^.  For  low  values  of  x,  fCx)^/  x  and  Eq. 

63  reduces  to 


— tn  T4  (u  .  Au  .)  ^ 
^  J  J  J  ” 


“■tn  T, 


h  =  1 


(64) 


If  all  the  Xj^’s  are  high,  f (x)  ^  (2x/tt)^  and  Eq.  63  reduces  to 


-  ^2 
-An  T .  (v  .  Av.)  -v  Z  (-An  t.) 
iJ>J“  h=l  ^ 


(65) 


Equations  64  and  65  are  valid  in  the  low  x  and  high  x  regions 
respectively  and  the  x  value  need  not  be  known  for  these 
calculations . 

The  parameters  (7/d)^  and  x^  are  written  in  Eq.  58  and 
Eq.  59  in  terms  of  the  zonal  band  model  parameters  which  are 
dependent  on  the  pressure  and  temperature  in  the  following  way . 

(s/d)j_  =  (s°/d)^p^  (66) 

and 

(7/d)^  (7°/d)h  Ph  +  pS  (8^) 
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Where  Pj^  and  are  the  pressure  of  the  absorbing  gas  the 
pressure  of  the  broadening  gas  respectively  in  the  isothermal 
zone  h  of  path  length  The  quantities  (s°/d)j^  and 

are  functions  of  temperature.  The  parameters  (7/d)^  and  may 
now  be  written  in  terms  of  partial  pressures  and  temperature 
dependent  band  model  parameters  as, 


1 

2 
i 


(68) 


and 


t 

[f  .  ^1>] 


(69) 


2^  (s°/d)^  p^“ 


Ph^  +(^'’b/‘»hPh 


For  the  case  of  several  broadening  gases  the  single  terms 

(7°^j/d)j^  Pj^^  for  a  given  zone  are  replaced  by  sums  over  species 

in  that  zone,  namely,  2  h  Ph^’ 

s 

For  the  spectral  region  where  two  gas  species  both  absorb 
radiant  energy,  the  transmittance  is  first  computed  for  each 
gas  separately  and  then  the  product  determines  the  effective 
transmittance.  For  example,  the  case  of  overlapping  H2O  and 
CO2  bands  near  2.85  [im  would  be  calculated  by 

HoO  C0„ 

?l(Vj,AVj)  -  (Vj.AVj)  (70) 

In  a  case  like  this  the  second  gas  is  always  treated  as  a 
broadening  agent  when  the  first  is  treated  as  the  absorber,  and 
vice  versa. 


31 


AEDC-TR-71-78 


The  Equation  of  Transfer 

The  infrared  brightness  method  of  obtaining  information 
about  the  temperature  distribution  in  a  hot  gas  from  its 
emission  and  absorption  spectra  was  suggested  by  Silverman 
(Ref.  107)  and  has  been  developed  and  verified  by  numerous 
investigators .  The  technique  is  based  on  the  recognition  that 
radiation  is  both  emitted  and  absorbed  in  a  hot  gas  in  local 
thermodynamic  equilibrium  and  that  the  spectral  irradiance  by 
a  one  dimensional  beam  of  monochromatic  radiation  varies 
according  to  the  equation  of  transfer. 

dH(Vj,s)  =  -k(Vj,s)  p(s)  [  H(vj,s)  -  N^j(Vj,Tg)  ]  ds.  (71) 
H(v.,s)  is  the  irradiance;  N,  (v.,T  )  is  the  Planck  function; 

J  D  J  S 

k(v.  is  the  spectral  absorption  coefficient;  p(s)  and  T 
J  f  s 

are  the  partial  pressure  and  temperature  of  the  gas  at  the 

frequency  Vj  and  the  position  s.  The  left  hand  side  of  this 

equation  is  the  change  in  Irradiance  on  going  from  a  point  s 
to  a  point  s  +  ds.  The  first  term  on  the  right  hand  side  is 
the  decrease  in  irradiance  due  to  absorption  along  the  path  ds; 
the  second  term  is  the  increase  in  irradiance  due  to  the 
emittance  of  the  hot  gas  within  ds.  The  formal  solution  of  the 
equation  of  transfer  may  be  written 

s 

H(v. ,s)  =  H(v^,0)  exp[-  /  kpds] 

J  J  Q 

+  J  N.  (v.,T  ,)exp[-  /®kpds]kpds  (72) 
o  s' 

where  the  optical  thickness  between  s  and  s'  is 

s 

t  =  f  k  p  ds  (73) 

s' 
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At  this  point  it  was  increasingly  evident  that  the  math¬ 
ematical  tools  and  experimental  techniques  were  available  for 
predicting  the  radiances  of  hot  gases,  even  with  temperature, 
pr^essure,  and  concentration  gradients,  provided  these  gas 
parameters  are  known  for  the  physical  situation  under  investi¬ 
gation  and  provided  sufficiently  accurate  data  on  the  spectral 
parameters  of  line  strength  and  line  width  are  available  to  use 
in  the  band  model.  It  was  also  quite  evident  that  very  little 
was  known  regarding  the  inverse  problem,  namely,  determining 
the  temperature,  pressure,  and  concentration  gradients  in  a  hot 
inhomogeneous  gas  from  spectroscopic  data  and  a  suitable  band 
model.  A  very  pertinent  effort  in  this  direction  was  that  of 
Dr.  William  F.  Herget  and  associates  (Ref.  52-Ref.  55)  who 
extended  the  spectroscopic  emission-absorption  method  to  obtain 
information  about  both  temperature  distribution  and  concentration 
distribution.  To  determine  these  quantities  the  spectral 
radiance  and  spectral  transmittance  of  the  infrared  active  spe¬ 
cies  in  the  hot  gas  must  be  measured  spectroscopically  along  a 
series  of  coplanar  chordal  lines  of  sight .  These  data  allow  the 
calculation  of  the  radial  distributions  of  temperature  and 
partial  pressure  of  the  species  of  interest  by  inversion  methods. 
The  method  is  best  explained  by  first  considering  media 
homogeneous  in  temperature  and  pressure,  and  then  extending  the 
technique  to  heterogeneous  media. 

The  partial  pressure  of  the  gaseous  medium  may  be  determined 
by  its  absorptive  properties  which  are  expressed  by  the  spectral 
transmittance  of  the  gas 

T(v)  = 

where  k(v,T) ,  the  spectral  absorption  coefficient,  is  a  unique 
property  of  the  absorbing  gas  for  a  given  frequency  and  temp¬ 
erature.  k(v,T)  can  be  completely  determined  for  a  particular 
species  by  performing  laboratory  absorption  measurements  over  a 
range  of  temperatures . 

The  spectroscopic  determination  of  partial  pressures  makes 
use  of  the  fact  that  once  k(v,T)  has  been  determined  for  various 
species  by  laboratory  measurements  then  the  partial  pressure  of  an 
active  gas  can  be  determined  from  Eq.  74  by  a  measurement  of  t(v), 
T  and  -t . 

The  temperature  may  be  determined  by  a  measurement  of  the 
spectral  radiance  which  is  the  power  of  the  radiation  per  unit 
frequency  interval  emitted  from  a  unit  area  of  the  gas  surface 
filling  a  given  solid  angle  and  is  expressed  in  units  of 
W[cm2  cm“^  str ]“■*•.  Spectral  radiance  is  measured  by  comparing 
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the  brightness  of  the  hot  gas  with  the  brightness  of  a  standard 
source  such  as  a  blackbody  at  a  known  temperature.  The  optical 
system  is  arranged  to  insure  that  the  detector  views  the  same 
solid  angle,  frequency  interval,  and  source  area  whether 
viewing  the  hot  gas  or  the  blackbody.  Because  of  this  arrange¬ 
ment  the  ratio  of  the  spectral  radiance  of  the  gas  to  the 
spectral  radiance  of  the  blackbody  is  the  same  as  the  ratio  of 
the  spectral  irradiance  of  the  gas  to  the  spectral  irradiance 
of  the  blackbody.  That  is 

N(v)/Njj(v)  =  H(v)/H^(v)  (75) 


Also  with  this  optical  arrangement  the  spectral  radiance  of  the 
hot  gas  will  be  equal  to  that  of  the  blackbody  if  the  detector 
produces  the  same  output  signal  no  matter  which  source  is  view¬ 
ed.  The  spectral  radiance  of  a  blackbody  Nb(v,T)  is  a  unique 
function  of  frequency  and  temperature  given  by  the  Planck  formula 


Nb(v,T)  =  2c^hv^[ exp(hcv/kT)-l]“^ 

=  1.1909  X  10" exp (1.439  v/T)-l]“^ 


(76) 


The  spectral  radiance’  of  a  homogeneous  gas  N(v)  is  related  to 
the  spectral  radiance  of  a  blackbody  at  the  same  temperature  by 


N(v)  «  Njj(v,T)  e(v) 

=  N^,(v,T)[1  -  T(v)]  (77) 


Thus  once  N(v)  and  t(v)  have  been  measured,  N^(v,T)  can  be 
calculated  from  Eq.  77  and  the  gas  temperature  is  then  uniquely 
determined  by  the  Planck  formula  since  the  temperature  is  then 
the  only  unknown  quantity  in  that  expression  and  is  equal  to  the 
gas  temperature. 

In  extending  measurements  to  an  unhomogeneous  gas  in  which 
temperature  and  pressure  gradients  exist  it  is  generally  assumed 
that  the  gas  consists  of  a  series  of  regions  each  of  which  is 
isothermal  and  isobaric  within  the  precision  of  measurement.  If 
the  gas  under  investigation  is  in  an  axially  symmetric  flow 
field  as  is  often  the  case  for  engine  exhausts  then  the  zones  of 
constant  temperature  and  pressure  are  chosen  as  shown  in  Figure 
4.  If  the  temperature  and  concentration  gradients  are  steep  the 
regions  would  have  to  be  small  but  in  principle  such  a  division 
can  always  be  made .  The  medium  between  the  sample  and  the 
detector  is  considered  zone  zero.  Spectroscopic  measurements 
are  made  along  lines  of  sight  as  shown  in  the  diagram. 

Herget  made  line  of  sight  measurements  of  medium  resolution 
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of  spectral  radiance  and  spectral  transmittance  at  small  fre¬ 
quency  intervals  centered  about  3  specific  frequencies  in  the 
infrared  for  gaseous  carbon  dioxide  in  a  rocket  exhaust.  The 
partial  pressures  and  temperatures  for  each  radial  zone  were 
determined  by  solving  appropriate  spectral  radiance  and  trans¬ 
mittance  equations  utilizing  an  "onion  peel"  or  back  substitu¬ 
tion  calculation.  This  means  that  the  partial  pressure  and 
temperature  for  the  outermost  zone  are  determined  directly 
from  line  of  sight  data  taken  through  that  zone  alone.  These 
results  are  then  used  with  the  data  taken  along  the  next  inner¬ 
most  line  of  sight  to  obtain  temperatures  and  pressures  for  the 
next  innermost  zone.  This  method  is  repeated  until  the  problem 
is  solved. 

Progress  in  developing  a  low  resolution  band  model  method 
for  determining  the  temperature,  pressure,  and  concentration 
gradients  in  a  hot  Inhomogeneous  gas  required  that  selections  be 
made  among  the  band  models,  the  available  approximations,  the 
spectral  region  and  the  experimental  techniques.  The  decision 
was  made  to  concentrate  on  the  work  of  two  groups  which  have 
been  very  active  in  the  development  and  application  of  band 
model  methods  to  predicting  radiative  heat  transfer  where  the 
concentrations  and  temperatures  of  combustion  gases  are  known. 
These  are  the  General  Dynamics  Convair  Group  of  C.  B.  Ludwig, 
J.A.L.  Thomson,  and  W.  Malkmus,  and  the  Warner  and  Swasey  Group 
of  H.  Babrov ,  R.  Tour in,  and  B.  Krakow.  The  data  of  the  General 
Dynamics  group  on  band  model  parameters  in  the  Infared  at  2.7^, 
were  to  provide  a  computional  base .  The  Warner  and  Swasey 
development  of  the  Curtis-Godson  approximation  utilizing  the 
statistical  band  model  with  Lorentz  lines  of  equal  width  and 
intensity  provided  the  mathematical  model.  This  band  model  and 
intensity  distribution  were  chosen  because  they  have  been  shown 
to  satisfactorily  and  conveniently  represent  the  water  vapor 
band  spectra  and  because  much  of  the  available  band  model  data 
were  based  on  the  statistical  model.  A  different  band  model, 
for  example  the  quasi-random  model,  would  possibly  provide  an 
even  better  representation  of  the  water  vapor  spectra  but  the 
computational  complexity  appeared  unnecessary  for  the  precision 
of  the  measurements  contemplated.  Although  several  approxi¬ 
mations  are  available  for  dealing  with  mixtures  of  gases  only 
the  Curtis-Godson  approximation  was  given  serious  consideration 
in  the  present  work  because  the  mathematical  development  of  this 
approximation  had  been  well  advanced  and  its  representation  of 
inhomogeneous  gases  under  operating  conditions  similar  to  those 
encountered  in  a  test  cell  environment  had  been  experimentally 
verified  by  the  Warner  and  Swasey  group. 

The  task  of  "inverting"  the  equations  of  the  Warner  and 
Swasey  group  was  undertaken  by  Mr.  LeRoy  Brewer,  ETF  Special 
Projects  Branch,  Arnold  Engineering  Development  Center.  The 
equations  referred  to  are  those  derived  in  this  report  in  the 
section  on  the  Curtis-Godson  approximation.  The  equations  were 
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originally  developed  to  predict  radiance  knowing  temperature 
and  concentration  profiles  of  inhomogeneous  gases.  The 
"inversion”  required  determining  temperature  and  concentration 
profiles  from  measured  values  of  spectral  radiance  and  spectral 
transmittance.  This  was  achieved  by  a  back  substitution 
computational  technique  utilizing  data  obtained  from  spectros¬ 
copic  emission-absorption  experiments.  The  band  model  develop¬ 
ment  by  Brewer  requires  the  tabulated  values  of  the  band  model 
parameters  So/d  and  7°/d  as  functions  of  temperature  compiled 
by  the  General  Dynamics  group  for  various  temperatures  ranging 
from  300  to  3000°K.  In  the  development  which  follows  the  geom¬ 
etry  of  the  problem  is  shown  in  Figure  4  which  was  discussed 
above  in  the  section  on  emission-absorption  measurements. 

INVERSION  OF  THE  CURTIS-GODSON  EQUATIONS  FOR  STATISTICAL  MODEL 

Given  the  following  expression  for  the  spectral  radiance 


N<i3.  An.)  = 
J  >  J 


n 
2 
i=l 


- 


T. (n .  An  )] 

J*  J  >  J 


(49) 


SECTOR  A. 


I.  In  sector  A,  this  equation  becomes  simply  the  equation  for 
Zone  1  and  thus, 


N,  (.v  Au.)  = 

J-  Of  J 


)tl  - 


AO.)] 

f  J 


where  the  transmittance  of  zone  zero  is  considered  to  be  unity. 

Once  N.  (n.  An.)  and  t.,  (n  .  An.)  have  been  measured,  N.  (n  .  T, ) 
Aj,j  J-JjJ  DJjJ. 

can  be  calculated  directly  from  this  equation  and  the  gas 
temperature  in  zone  1  is  then  uniquely  determined  by  the  Planck 
function 


NbCnfT)  =  2c^hn^  [  exp  (hen AT) -1]"^  (76) 

Result :  T J 

II.  The  band  model  parameters  (S°/d).  and  (7°/d).  are  tabulated 
fuctlons  of  temperature  and  can  now  be  determined  for  zone  1. 

Result:  (S°/d)j^  and  (7°/d)j^ 

t  ]k 

III.  Given  the  definition  of  x.  for  the  h  isothermal  zone 
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(S/d)j^^h  (S°/d)jj 

^  2t(7  /d)jj  2ir(7°/d)uPu 

The  pressure  term  cancels;  For  zone  1  the  pathlength  is 
measured  and  the  values  of  (SO/d)^^  and  (7°/d)j^  are  known. 

Result ; 

IV.  Given:  f (x) ,  the  Ladenburg  and  Reiche  function  the  values 
of  which  are  tabulated. 

Result:  f(x^) 

V.  Given:  The  transmittance  of  a  homogeneous  gas  for  a 
statistical  band  model  with  constant  line  widths 

-tn  T  ==  2ir(7/d)  f  (x) 

Then  for  Sector  A  which  is  zone  1 

-^n  =  2ir(7°/d)j^  Pj^f(Xj^) 

Here  everything  is  known  except  P,  for  which  a  direct 
solution  is  obtained 

Result : 

For  zone  1  we  now  know  the  temperature  and  the  partial 
pressure  of  the  absorbing  gas 

Sector  B 

I.  Measure  N^C-o.  t  (u  AiJ.)  =  t  „  (u  AD.) 

J I  J  "  jf  j  ji  j 


and  the  path  length  for  the  various  zones 
The  equation  for  the  spectral  radiance  becomes 


+  Nb(Vj_Ti)[?i_2(Vj_AVj)  - 
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There  are  two  unknowns  in  this 


T 


1,2 


equation  N.  (u.  T„)  and 

D  J  »  ^ 


II.  The  transmittance  of  the  inhomogeneous  path  in  terms  of  the 

/ 

zonal  parameters  is  given  by  Eq.  60  and  for  the  present 
case  takes  the  form, 


The  unknowns  in  this  equation  are  (-^,n  Tg)  ,  Xg,  and  fCXg) 
where  fCXg)  is  determined  from  Xg  by  the  Ladenburg  and 
Reiche  function. 

III.  These  two  equations  for  sector  B  contain  the  information 
for  determining  the  pressure  and  temperature  of  zone  2. 
However  the  equations  cannot  be  written  in  terms  of  these 
unknowns,  Pg  and  Tg,  and  solved  simultaneously  because 
of  the  Ladenburg  and  Reiche  function  and  because  S°/d 
and  7°/d  are  tabular  functions  of  T.  Therefore  as  an 
initial  step  in  an  iterative  solution  assume  one  of  the 
approximations,  either  the  low  x  or  the  high  x  approximations 
to  equation  60  is  valid.  For  example,  assuihe  the  validity 
of  the  low  X  or  Beer*s  law  approximation;  then, 

-^n  Tg  =  -In  ^12  1^  *^1  -  -tn  Tg  -  In 


38 


AEDC-TR-71-78 


and 

-In  2  “  '^1  ~  "^2 

The  first  of  these  equations  can  be  solved  for  -^n  Tg  since 
--tn  T,  „  T  is  measured  and  — tn  t,  is  known  from  sector  A. 

1  I  A  }  1  J-  ^ 

The  second  of  these  equations  then  yields  (-^n  t  g)  where 
the  prime  indicates  an  approximate  value  for  the  trans¬ 
mittance  through  zones  1  and  2 . 

IV.  Substituting  -^n  g  equation  78  is  solved  for  N^(^j^Tg') 
and  the  temperature  of  zone  2,  Tg',  is  obtained  directly 
from  the  Planck  radiation  expression. 

V.  Following  the  routine  of  sector  A,  the  band  model  parameters 
(S°/d)^  and  (7°/d)^  are  obtained  for  the  temperature  T^ 

and  lead  to  a  value  for  x^  through  Eq.  56,  and  thence 
to  fCx'g)  by  the  tabulated  Ladenburg  and  Reiche  function. 

VI.  Now  substitute  (-- tn  x'g,  along  with  the  known  zone 

1  parameters  and  calculate  (--tn  T'g).  There  are  two  possible 
consequences  of  this  procedure; 

a)  If  T'g  =  Tg,  then  Beer’s  law  is  valid  and  we  move  on  to 
Step  VII. 

b)  If  T’g  ^  Tg,  then  adjustments  must  be  made  in  the  values 
of  (-^n  'T'g),  x’g,  and  f(x’g).  The  adjustment  in  these 
parameters  may  be  accomplished  by  examining  the  other 
approximation  to  Eq.  60,  for  example,  the  high  x  approxi¬ 
mation  . 

VII.  Agreement  between  the  calculated  and  measured  transmittance 
of  Sector  B  is  taken  to  imply  the  correctness  of  the 
parameters  of  zone  2,  in  particular  the  temperature  Tg 
and  the  transmittance  Tg.  The  correct  temperature  yields 
correct  (S°/d)g  and  (7°/d)g  from  tables.  These  in  turn 
provide  correct  Xg  and  fCxg).  The  partial  pressure  of  the 
infrared  active  gas  in  zone  2  is  now  determined  from 
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--tn  Tg  =  2T(7°/d)2  Pg  fCXg) 

Sector  C. 

I.  Measure  N(,(»j,AUj)  .  “  '"l, 2, 3, 2,1 

and  the  path  lengths  for  the  various  zones.  The  equation 
for  the- spectral  radiance  becomes 

+  Nb(Vj,T2)['r.i(nj,APj)  -  **^1 , 2  j  » j  ^  ^ 

+  Nb(Vj,T3)[?i_2('=j.AUj)  -  ?i,2,3('5j,Aej): 

+  -  ’'1,2, 3, 2 

(Uj.AUj)] 

^b^^J  ^1^1^  ^1  2  3  "'*^1  2  3  2  1 

(ej.AUj)] 

There  are  three  unknowns  in  this  equation: 

Nb(Oj.T3).  ?a,2,3("j.^'=J>>  ““  ^l,2,3,2(”j,^'=j^ 

II.  The  transmittance  of  the  Inhomogeneous  path  in  terms  of  the 
zonal  parameters  takes  the  form  shown  on  the  following  page. 
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—tn 


/  “-tn 

V  /-^n  Tg> 

1  +  O-TT  1 

\  ’■3') 

1" 

I'^TxJT', 

/  +  ^''2Vf(Xg)  ; 

/  ^3\TTxJ7“; 

f 

— tn  T,  -tn  Tg  — tn 

^*1  f (Xj)'  +  ^^2  f(x^i  +  ^"3  f(x-y 


2x, 


— tn  r. 


-^n  T 


+  X, 


f  (x. 


?) 


2X-, 


V  f(X^)/  K  fCXg)/  HfCXg)  / 


There  are  three  unknowns  in  this  equation: 
(--tn  Tg),  Xg,  and  f(Xg) 


III.  As  before  one  assumes  an  approximation  to  Eq,  60  as  an 
initial  step  in  an  iterative  solution  for  the  pressure 
and  temperature  in  zone  3.  Assuming  Beer’s  law,  for 
example,  gives, 

-In  2, 3,2,1 


=  —tn  -  -tn  Tg  -  -tn  Tg  -  -tn  Tg  -  -tn 

and 

-In  ^12^3  2  ^  ”^1  -  T2  “  Tg 

“-tn  g  g  =  “-tn  T  “  -tn  Tg  “  -tn  Tg 


IV. 


The  solution  of  these  equations  leads  to  approximate  values, 

-tn  ^'123  ^12  3  2 

Substitute  “-tn  2  3  ^1  2  3  2 

solve  for  Nj^  ,  T'g)  .  The  approximation  to  the  temp“ 
erature  of  zone  3,  T'g,  is  obtained  from  the  Planck 
radiation  expression.  Following  the  routine  of  sectors  A 
8e  B  yields  an  approximate  value  for  x'g  and  for  the 
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Ladenburg  and  Reiche  function  fCx'^).  Substituting 
(-■ tn  T'g),  x'g  ,  and  fCx'g)  together  with  previously 
determined  quantities  into  £q.  80  yields  the  approximate 
transmittance  for  Sector  C, 

V.  As  before,  if  the  calculated  transmittance  for  Sector  C, 

T^,  equals  the  measured  transmittance  for  Sector  C,  t  , 
c  c 

then  the  Beer*s  law  approximation  is  valid.  If  however, 
x'  X  then  adjustment  in  the  parameters  is  necessary. 
When  sufficiently  precise  agreement  is  reached  between 
calculated  and  measured  transmittance  of  Sector  C  the 
parameters  of  zone  3  may  be  assumed  correct  and  the 
temperature  and  partial  pressure  of  the  gas  may  be 
determined  for  that  zone. 

Other  Sectors. 

The  complete  problem  is  solved  by  repeating  the  above 
procedure  for  each  sector  until  all  chosen  sectors  have  been 
evaluated. 

Presence  of  Foreign  Gases 

In  the  inversion  process  described  above  no  provision  was 
made  for  the  presence  of  any  gas  other  than  the  one  whose 
radiance  and  transmittance  is  being  observed.  In  most  com¬ 
bustion  processes  of  interest  this  situation  is  not  obtained. 
There  are  usually  broadening  agents  present  in  appreciable 
concentrations  and  the  foreign-gas  broadening  effects  and  band 
overlap  effects  must  be  considered.  The  band  model  halfwldth 
which  enters  into  the  radiance  equations  for  a  given  isothermal, 
isobaric ,  zone  within  the  inhomogeneous  gas  can  be  written 

(7/d)^  “  Ph*  +  (7b°/d)^  (81) 

where  the  index  a  indicates  the  absorbing-emitting  gas  and  the 
index  b  indicates  the  non-absorbing  or  broadening  gas; 

Pjj®’  and  pj^^  are  the  partial  pressures  of  the  absorbing  and 
broadening  gases  respectively  in  zone  h.  The  band  model  half- 
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Width  for  the  foreign  gas  can  be  approximated  in  terms  the  band 
model  halfwidth  of  the  absorbing-emitting  gas  by  writing 

=  K  C-VaVd)  (82) 

where  K  is  a  constant  which  differs  for  each  broadening  gas 
but  which  is  assumed  in  this  approximation  to  be  Independent 
of  temperature  and  frequency.  Also  observe  that  if  only  two 
gases  are  present  the  partial  pressure  of  the  broadening  gas 
may  be  expressed  in  terms  of  the  absorbing-emitting  gas  as 

p^  =  1  -  p^  (83) 

The  band  model  half width  parameter  then  takes  the  form, 

(7/d)^  .  (7^°/d)^[K  +  (1  -  K)  (84) 

For  the  case  of  several  broadening  gases  the  tem  (VK^/d)  j^Pi, " 

is  replaced  by  Z  (7^  ^h  summation  over  the 

s 

index  s  indicates  that  one  should  s\im  the  product  of  partial 
pressure  and  broadening  parameter  for  each  species  for  each 
Zone  h.  Although  this  additional  complication  can  be  incor¬ 
porated  into  the  equations,  it  is  in  most  cases  not  necessary 
to  consider  the  effect  of  more  than  one  broadening  gas.  The 
case  of  interest  here  is  that  in  which  the  absorbing-emitting 
gas  is  water  vapor  and  the  broadening  gas  is  carbon  dioxide 
or  perhaps  nitrogen.  There  a  few  experimental  data  from  which 
a  value  of  K  may  be  obtained.  Babrov  (Ref.  73)  estimated  that 
for  water  vapor  lines  broadened  by  nitrogen  ICj  is  approximately 

0.2.  This  is  in  agreement  with  the  General  Dynamics  Group  (Ref. 
44)  who  give  a  value  7^  =:  0.09  when  the  water  vapor  parameter 

is  7tt  n  0.44  for  nitrogen  broadening  water  vapor  lines.  In 

exhaust  gases  CO2  is  the  primary  foreign  gas  in  H2O  absorption 
studies  and  has  a  much  greater  broadening  ability  than  N2. 
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Rusk  (Ref.  106)  found  CO2  to  be  1.6  times  as  effective  as 
nitrogen  in  broadening  water  vapor  lines.  The  data  of  the 
General  Dynamics  group  support  this  value  also.  On  the 
basis  of  these  findings  it  is  suggested  that  an  initial  value 
of  0.32  be  assigned  to 

Very  little  experimental  information  is  actually  available 
regarding  half widths.  This  is  true  not  only  for  band  models, 
but  for  the  basic  spectral  parameter  single  line  halfwidths 
for  pure  homogeneous  gases  as  well.  The  self  broadening  and 
foreign  gas  broadening  parameters  as  functions  of  frequency 
and  temperature  are  not  known  for  most  molecules  and  should 
be  measured  for  those  molecules  of  importance  in  combustion. 

The  spectral  region  where  the  measurements  of  interest 

to  this  report  are  made,  namely,  the  2.7|j,  infrared  region  is  a 

region  where  both  water  vapor  and  carbon  dioxide  are  absorbing 

and  emitting  radiation.  When  this  is  the  case  for  two  gas 

species,  more  accurate  values  of  the  spectral  radiance  and  the 

spectral  transmittance  are  obtained  if  one  first  computes 

T.(n.  At).)  for  each  absorbing  gas  separately,  and  then  multl- 
1  J »  3 

plies  the  two  t.  (u  At).)  s  to  obtain  the  effective  transmittance. 

1  J  >  J 

In  this  case  for  example, 

T.(Oj  Ai3j)  =  [T.®2®][i^^°2]  (85) 

In  these  calculations  the  second  gas  is  treated  as  the  broadening 
gas  when  the  first  is  the  absorbing  gas  and  vice  versa.  This 
refinement  can  be  incorporated  into  the  calculations  for  temp¬ 
eratures  and  partial  pressures,  but  the  additional  complication 
is  not  warranted  at  this  time  by  the  precision  of  the  experiments. 

Another  point  which  should  be  considered  with  care  in  the 
development  and  utilization  of  the  inversion  process  is  the 
transfer  of  the  zonal  spectral  transmittance,  j  * j ^  ’  from 

one  sector  to  another.  The  optical  depth  for  a  given  zone 
differs  from  sector  to  sector  as  is  seen  by  examining  Fig.  4. 
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It  is  possible  that  the  path  length  in  zone  1  sector  A  may  be 
greater  than  the  combined  path  lengths  through  both  pieces  of 
zone  1  in  sector  B.  A  similar  situation  exists  for  the  other 
zones  although  it  becomes  less  pronounced  as  the  center  zone 
is  approached.  Perhaps  the  simplist  way  to  account  for  this 
situation  is  to  observe  that  the  logarithm  of  the  spectral 
transmittance  may  be  written 

-  -tn  T  =  k(^^T)  p-t  (86) 

where  t  is  the  optical  path  length  in  the  medium.  Then  if  ^ 
is  reduced  to  one-half,  the  t  is  also  reduced  to  one-half, 
provided  pressure  is  not  changed  appreciably. 

Calculations 

An  examination  of  the  utility  of  the  inversion  procedure 
was  initiated  using  a  desk  calculator.  Calculations  were  made 
for  a  two  zone  situation  in  which  the  temperature  and  partial 
pressure  differ  for  the  two  zones.  The  data  input  for  the 
calculation  were  taken  by  the  Warner  and  Swasey  group  who 
measured  transmittance  as  a  function  of  temperature  and  pressure 
using  furnaces  under  controlled  laboratory  conditions.  Applying 
these  data  to  the  Inversion  procedure  utilizing  a  statistical 
band  model  with  the  Curtis-Godson  approximation  the  temperature 
and  pressure  of  each  zone  are  obtained.  For  zone  one  the 
calculated  pressure  is  within  2%  of  the  measured  pressured  and 
the  calculated  temperature  within  1%  of  the  measured  temperature. 
For  zone  2,  the  agreement  was  within  4%  for  pressure,  1%  for 
temperature,  and  3%  for  transmittance.  This  indicates  self 
consistency  of  the  Inversion  process.  When  the  band  model 
parameters  S^/d  and  as  functions  of  temperature  are 

taken  from  tables  compiled  by  the  General  Dynamics  group. 

The  error  in  pressure  calculations  is  slightly  higher  than 
50%.  This  discrepancy  was  traced  to  a  difference  in  the  value 
of  7°/d  determined  by  averaging  and  by  direct  measurement  over 
a  small  wavelength  interval. 

Two  other  points  of  interest  were  brought  out  by  the  hand 
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calculations.  Because  the  Beer-Lambert  lav  Is  not  valid,  the 
calculation  requires  an  iteration.  The  calculation  verified 
that  this  iteration  process  converges  very  rapidly  to 
acceptable  values  of  pressure  and  temperature.  The  calculation 
also  revealed  that  the  high  x  approximation,  namely, 

-^n  T  [2  (-^n 

^  h  ^ 

provides  more  correct  values  to  the  transmittance  and  thus 
more  rapid  convergence  than  does  the  low  x  approximation  for 
gas  samples  in  the  temperature  and  pressure  range  of  Interest . 
(600  to  1200“K  and  100  torr) 
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SECTION  V 
EXPERIMENTAL 


Two  types  of  experiments  were  designed  in  cooperation  with 
the  Special  Projects  Group,  ETF,  Arnold  Engineering  Development 
Center  for  obtaining  data  under  controlled  laboratory  conditions. 
The  first  of  these  had  the  objective  of  assessing  the  utility 
and  validity  of  the  previously  formulated  band  model  inversion 
technique  as  a  diagnostic  tool  for  determining  temperatures  and 
partial  pressures  from  spectroscopic  measurements  of  radiance 
and  transmittance.  The  second  experiment  was  designed  for 
measuring  the  effect  of  foreign  gas  broadening  on  the  band  model 
half width  parameter. 

Only  a  few  comments  will  be  made  regarding  the  second 
experiment  because  there  was  not  time  to  carry  it  out  under  the 
auspices  of  the  present  contract.  This  would  be  an  absorption 
experiment  in  which  band-width  measurements  are  made  as  a 
function  of  temperature  and  partial  pressure  on  hot  water  vapor 
mixed  with  such  gases  as  carbon  monoxide,  carbon  dioxide,  or 
nitrogen.  These  measurements  can  be  accomplished  satisfactorily 
with  an  instrument  of  low  resolution  such  as  a  Perkin-Elmer 
Model  12  infrared  spectrometer  fitted  with  a  prism  or  grating. 

The  difficulties  of  such  measurements  are  not  optical  but  are 
associated  with  the  gas  handling  system  and  with  the  removal  of 
atmospheric  water  vapor  from  the  optical  path.  The  later 
difficulty  is  most  satisfactorily  solved  by  placing  the  entire 
experiment  in  a  vacuum.  This  is  probably  not  possible  in  the 
existing  test  cell  environment  and  so  the  purging  of  the  system 
with  an  inert  gas  would  be  required.  The  problems  associated 
with  the  gas  handling  system  for  such  an  absorption  experiment 
are  more  restrictive.  An  absorption  cell  is  required  which 
remains  vacuum  tight  as  temperatures  range  from  300®K  to  1200® K 
and  back  to  300®K.  The  overall  gas  handling  system  must  provide 
for  the  mixture  of  two  or  more  gases  at  elevated  temperatures, 
the  accurate  measurement  of  their  partial  pressures,  and  the 
transfer  of  the  mixture  to  the  absorption  cell.  Because  of 
these  problems  consideration  is  .being  given  to  design  of  a  com¬ 
bustion  experiment  for  measuring  foreign  gas  broadening  of  the 
molecular  band. 

The  objective  of  the  first  experiment  was  of  more  immediate 
concern  and  the  decision  was  made  to  seek  emission  and  absorption 
data  at  2.7|i  in  the  spectrum  of  hot  water  vapor  produced  in 
flames.  These  experiments  were  carried  out  by  UTSI  personnel  in 
cooperation  with  the  ETF  Special  Projects  Group  utilizing  the 
special  burners  and  spectrometers  available  at  AEDC. 

To  provide  the  required  llne-of-sight  radiance  and  trans¬ 
mittance  data  the  infrared  instrumentation  system  must  be 
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capable  of  the  following  experiments; 

1.  Measurement  of  flame  spectral  radiance  is 
made  by  comparing  flame  intensity  with 
blackbody  intensity;  radiation  is  optically 
chopped  between  the  flame  and  the  detector. 

2.  Measurement  of  flame  spectral  transmittance 
is  made  by  observing  the  radiation  from  a 
greybody  source  after  passage  through  the 
flame;  greybody  radiation  is  optically 
chopped  between  the  greybody  and  the  flame. 

3 .  The  line  of  sight  miist  be  varied  by  spatially 
scanning  the  flame  image  across  the  entrance 
slit  of  the  spectrometer. 

Figure  5  shows  the  emission-absorption  experimental 
arrangement . 

The  experiment  was  set  up  in  the  RTF  test  area  of  the 
Arnold  Engineering  Development  Center.  A  burner  whose 
characteristics  were  well  known  from  earlier  experiments  was 
made  available.  The  selection  of  a  burner  with  known 
characteristics  was  necessary  to  provide  a  comparison  standard 
for  the  spectroscopic  measurements.  Temperature,  partial 
pressure,  and  velocity  profiles  had  been  mapped  for  this  burner 
for  several  different  fuel-oxidizer  ratios  and  inlet  pressures. 
The  burner  was  placed  in  position  so  that  a  vertical  flame  would 
be  obtained.  Traversing  mechanisms  were  Installed  so  that  the 
flame  could  be  moved  across  the  optical  path  in  both  a  horizon¬ 
tal  and  a  vertical  plane.  Measurements  can  consequently  be  made 
at  various  stations  downstream  of  the  flameholder  and  trans¬ 
versely  at  each  station.  A  hood  and  exhaust  system  were  des¬ 
igned  and  placed  in  position  above  the  flame.  This  is  necessary 
because  the  major  portion  of  the  exhaust  gas  will  be  water  vapor 
and  this  must  be  removed  from  the  vicinity  of  the  spectroscopic 
measurements.  The  gas  handling  apparatus  was  designed  for  burn¬ 
ing  hydrogen  premixed  in  either  air  or  oxygen.  The  initial 
measurements  were  made  of  hydrogen  burning  in  air.  The  flame  was 
surrounded  by  a  shroud  gas  of  nitrogen.  On  occasion  argon  or 
heliiun  was  used  as  the  shroud  gas. 

A  Perkin-Elmer  model  112  double  pass  prism  spectrometer  for 
the  infrared  was  acquired  and  positioned  for  obtaining  spectral 
data.  This  instrument  utilizes  a  lead  sulfide  detector,  a  I3cps 
chopper,  a  preamplifier,  and  a  Perkin-Elmer  model  107  amplifier 
for  receiving  emission  radiation  directly  from  the  flame.  Across 
the  flame  directly  opposite  the  spectrometer  a  graybody  radiation 
source  has  been  placed  for  making  the  absorption  measurements. 
This  radiation  source  is  a  water  cooled  Globar  which  has  high 
emissivity  throughout  the  near  infrared  and  is  operated  at  about 
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1200° K.  A  movable  shutter  between  the  Globar  and  the  flame 
allows  the  Globar  to  be  isolated  so  that  Its  radiation  is  not 
received  at  the  spectrometer  when  it  is  desired  to  look  only  at 
the  flame.  The  signal  detection  system  for  the  Globar  consists 
of  a  300  cps  chopper  together  with  a  lock-in  amplifier.  The 
entrance  slit  of  the  spectrometer  is  about  18  Inches  from  the 
flame  on  one  side  and  the  Globar  is  the  same  distance  from  the 
flame  on  the  other  side.  An  optical  system  external  to  the 
spectrometer  connects  the  two,  placing  a  one-to-one  image  of 
the  Globar  on  the  center  line  of  the  flame  and  then  a  one-to-one 
image  of  the  Globar  image  onto  the  entrance  slit . 

The  optical  paths  between  the  Globar  and  the  flame  and 
between  the  spectrometer  and  the  flame  are  enclosed  in  trans¬ 
parent  plastic  boxes.  When  measurements  are  being  made,  the 
spectrometer  and  much  of  the  optical  path  is  flushed  with  an 
inert  gas  to  reduce  absorption  by  atmospheric  water  vapor .  This 
flushing  gas  will  generally  be  the  same  as  the  shroud  gas,  that 
is,  nitrogen,  argon,  or  helium. 

Preliminary  measurements  of  the  emission  spectra  of  the 
hydrogen-air  flame  were  made  in  the  spectral  region  1.5p,  to  2.8p,, 
showing  a  strong  molecular  band  at  approximately  2,5[i  and  a 
weaker  band  in  the  neighborhood  of  1.8p,.  The  structure  at  2,5[i 
consists  of  the  water  vapor  band  of  interest  and  the  first  over¬ 
tone  band  of  C02<  The  band  at  1.8[x  is  also  water  vapor. 

Preliminary  absorption  measurements  were  also  made  on  the 
flame  center  line  throughout  this  spectral  region  using  the 
Globar  as  a  greybody  source.  These  measurements  showed  a  very 
weak  absorption  in  the  region  of  the  water  vapor  band  at  2.5p,. 

Following  these  preliminary  emission-absorption  measure¬ 
ments  the  instrument  was  calibrated  using  a  secondary  standard 
tungsten  ribbon  source.  This  together  with  the  emission- 
absorption  data  gave  an  estimate  of  the  flame  temperature  which 
was  within  20%  of  the  previously  determined  probe  values  under 
similar  burning  conditions. 

Emission-absorption  measurements  were  then  performed  at  5mm 
intervals  across  the  flame  to  obtain  more  detailed  information. 
These  data  provide  the  input  to  the  band  model  inversion  program 
from  which  temperature  profiles  and  water  vapor  partial  pressure 
profiles  across  the  flame  could  be  calculated. 

Table  1  through  Table  IV  present  typical  sets  of  data  for 
four  different  stations  downstream  of  the  burner  base.  At  each 
station  the  flame  was  scanned  transversely  from  edge  to  edge  at 
5  mm  intervals  resulting  in  twenty-one  measurements  of  radiance 
and  transmittance  at  each  station.  These  data  were  then  invert- 
.ed  to  provide  temperatures  and  water  vapor  partial  pressures  for 
each  of  ten  zones  for  every  station  downstream. 
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Figure  6  presents  typical  spectral  radiance  and  trans¬ 
missivity  curves  across  the  flame  from  centerline  to  edge. 

In  Figure  7  typical  profiles  are  shown  of  temperature  and 
water  vapor  partial  pressure.  The  temperature  and  pressure 
profiles  obtained  from  a  chemical  equilibrium  calculation 
and  previous  probe  measurements  are  also  shown  in  Figure  7. 
This  profile  is  distinguished  from  those  determined 
spectroscopically  in  two  respects .  First ,  it  is  flat 
showing  an  cylindrical  isothermal  region  of  3/8  inch  radius 
about  the  centerline  of  the  flame.  Second,  its  value 
differs  from  the  spectroscopic  value  in  the  zone  of  maximum 
value  by  approximately  12%.  It  should  be  observed,  however, 
that  the  spectroscopic  values  are  fairly  high  but  show 
gradients  in  the  region  shown  by  the  probe  to  be  isothermal. 

The  reason  for  these  discrepancies  is  not  known.  The 
fault  may  rest  with  the  experiment  or  in  the  computer  program 
and  efforts  are  continuing  to  determine  which. 
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SECTION  VI 

CONCLUDING  REMARKS 

The  investigation  of  the  feasibility  of  applying  the  theory 
of  band  models  to  the  diagnostic  problem  of  determining  temp¬ 
eratures  and  concentrations  has  lead  to  the  development  of  a 
mathematical  inversion  technique  based  on  low -resolution  emission 
absorption  spectroscopy.  Band  model  theory  and  much  of  its 
literature  is  summarized  in  this  report .  The  inversion  of  the 
Curtls-Godson  approximation  of  a  statistical  band  model  with 
Lorentz  lines  of  constant  width  is  presented  in  detail.  Exper¬ 
iments  were  performed  to  verify  the  technique  but  are  incon¬ 
clusive  at  this  time. 

The  present  work  is  not  complete  and  its  goal  will  be 
reached  only  if  it  stimulates  further  effort  in  the  application 
of  spectroscopic  band  models  to  the  diagnostic  problem.  The 
present  band  model  inversion  needs  to  be  refined  and  the  utility 
of  other  band  models  should  be  investigated.  The  application  to 
molecules. other  than  H2O,  in  particular  COnt  and  to  other 
spectral  regions,  such^as  the  15|j.  region,  should  be  pushed. 

Both  foreign  gas  broadening  and  self  broadening  are  too  little 
understood  and  available  data  all  too  meager.  The  band  model 
parameters  are  not  known  for  enough  molecules,  over  an 
adequate  spectral  range,  with  sufficient  accuracy  for  really 
satisfactory  solution  of  the  diagnostic  problem.  Much 
fundamental  work  needs  to  be  done  in  these  areas.  It  is  not 
unreasonable  to  expect  the  final  goal  of  band  model  diagnostics 
to  be  on-line'  instrumentation  functioning  automatically  to 
provide  instantaneous  temperature  and  concentration  profiles 
for  the  gases  of  interest . 


51 


AEDC-TR-71-78 


BIBLIOGRAPHY 

1.  Elsasser,  W.  M.,  Phys.  Rev.  54,  126  (1938)  The  Absorption 
for  a  Band  Composed  of  IdentTcal  Spectral  Lines.  Uniformly 
Spaced . 

2.  Elsasser,  W.  M. ,  "Heat  Transfer  by  Infrared  Radiation  in 
the  Atmosphere",  Harvard  Meteorological  Studies  No.  6., 
Harvard  University  Press,  (1942). 

3.  Goody,  R.  M.,  Quart.  J.  Roy.  Meteorol.  Soc .  78,  165 
(1952)  Statistical  Band  Model. 

4.  Goody,  R.  M. ,  "Atmospheric  Radiation"  Oxford  University 
Press,  New  York,  (1964). 

5.  Wolfe,  \W.  L.,  Editor,  "Handbook  of  Military  Infrared 
Technology",  Office  of  Naval  Research,  Department  of  the 
Navy  (1965)  for  sale  by  Superintendent  of  Documents, 

U.  S.  Government  Printing  Office,  Washington,  D.  C. 

6.  Anding,  D.,  NAVSO  P-2499-1  (Feb.  1967),  AD-815481, 

Band  Model  Methods  for  Computing  Atmospheric  Slant  -  Path 
Molecular  Absorption. 

7.  Rochelle,  W.  C.,  NASA  TMX-53579  (Feb.  1967),  Review  of 
Thermal  Radiation  from  Liquid  and  Solid  Propellant  Rocket 
Exhausts . 

8.  Benedict,  W.  S.,  Herman,  R.,  Moore,  G.  E.,  and  Silverman, 
S.,  Can.  J.  Phys.  830,  (1956),  The  Strengths,  Widths, 
and  Shapes  of  Infrared  Lines.  1.  General  Considerations. 

9.  Benedict,  W.  S.,  Herman,  R.,  Moore,  G.  E.,  and  Silverman, 
S.,  Can.  J.  Phys.  850  (1956).  The  Strengths,  Widths, 
and  Shapes  of  Infrared  Lines.  2.  The  HCI  Fundamental. 

10.  Mason,  A.  A.,  and  Nielsen,  A.  H.,  J.  Opt.  Soc.  Am.  57 , 

1464  (1967)  Rotational  Spectrum  of  Hydrogen  Fluoride: 
Frequencies  and  Linewidths . 

11.  Curtis,  A.  R.,  Quart.  J.  Roy.  Meteoral.  Soc.  78,  638 

(1952)  .  — 

12.  Godson,  W.  L.,  J.  Meteor  ]^,  272  (1955). 

13.  Plass,  G.  N.,  J.  Opt.  Soc.  Am.  42,  677  (1952).  A  Method 
for  the  Determination  of  AtmospEeric  Transmission  Functions 
from  Laboratory  Absorption  Measurements  Applied  to  Atmos¬ 
pheric  Heat  Transfer  Problem. 

14.  Plass,  G.  N.,  J.  Opt.  Soc.  Am.,  690  (1958)  Models 
for  Spectral  Band  Absorption. 


52 


AEDC-TR-71-78 


15. 

16. 

17. 

18. 

19. 

20. 

21. 


22. 

23. 

24. 

25. 

26. 

27. 

28. 

29. 


Plass,  G.  N.  ,  J.  Opt.  Soc.  Am.,  821  (1959)-,  Spectral 

Emissivity  of  Carbon  Dioxide. 

Plass,  G.  N.,  J.  Opt.  Soc.  Am.  868  (1960)  Useful 

Representations  for  Measurements  of  Spectral  Band 
Absorption. 

Wyatt,  P.  J.,  Stull,  V.  R.,  and  Plass,  G.  N.,  J.  Opt.  Soc. 
Am.  1209  (1962),  Quasi-Random  Model  of  Band  Absorption. 

Plass,  G.  N.,  Appl.  Opt.  2,  515  (1963). 

Plass,  G.  N.,  Appl.  Opt.  4,  69  (1965).  Spectral  Band 
Absorptance  of  Radiation  Traversing  Two  or  More  Cells 
in  Series. 


Plass,  G.  N.,  Appl.  Opt.  4,  161  (1965).  The  Theory  of 
Absorption  of  Flame  Radiation  by  Molecular  Bands. 


Howard,  J.,  Burch,  D.,  and  Williams,  D.,  J.  Opt.  Soc. 
Am.  46,  (1956)  Infrared  Transmission  of  Synthetic 
Atmos^eres 


-186,  I 
-237,  II 
-242,  III 
-334,  IV 
-452,  V 


Instrumentation 
Absorption  by  Carbon  Dioxide 
Absorption  by  Walter  Vapor 
Application  of  Theoretical  Band  Models 
Absorptj,on  Laws  for  Overlapping  Bands. 


Burch,  D.  E.,  Singleton,  E.  B.,  and  Williams,  D.,  Appl. 
Opt.  1,  359  (1962). 


Burch,  D.  E.,  and  Williams,  D.,  Appl.  Opt.  1^,  473  (1962), 
Total  Absorptance  by  Nitrous  Oxide  Bands  in  the  Infrared. 


Burch,  D.  E.,  and  Williams,  D.,  Appl.  Opt.  1,  587  (1962), 
Total  Absorptance  of  Carbon  Monoxide  and  MeThane  in  the 
Infrared . 


Burch,  D.  E.,  Gryvnak,  D.,  and  Williams,  D.,  Appl. 
Opt.,  1^,  759  (1962),  Total  Absorptance  of  Carbon 
Dioxide  in  the  Infrared. 

Burch,  D.  E.,  and  Gryvnak,  D.  A.,  Appl.  Opt.,  2,  585, 
(1963).  “ 


Burch,  D.  E.,  France,  W.  L.,  and  Williams,  D.,  Appl. 

Opt.  2,  585  (1963). 

Burch,  D.  E.,  and  Williams,  D.,  Appl.  Optics  55  (1964), 
Tests  of  Theoretical  Absorption  Band  Model  Approximations. 

Burch,  D.  E.,  and  Gryvnak,  D.  A.,  JQSRT  6,  229  (1966). 
Laboratory  Investigation  of  the  Absorption  and  Emission 
of  Infrared  Radiation. 


53 


AEDC-TR-71-78 


30.  Burch,  D.  E.,  Gryvnak,  D.  A.,  and  Patty,  R.  R.,  J.  Opt. 

Soc .  Am.  57,  885  (1967).  Absorption  of  Infrared  Radiation 
by  COg  anHTlgO  Experimental  Techniques . 

31.  Burch,  D.  E.,  Gryvnak,  D.  A.,  and  Patty,  R.  R.,  J.  Opt. 

Soc.  Am.  58,  335  (1968).  Absorption  of  Infrared  Radiation 
by  COg  between  8000  and  10,000  cm“l. 

32.  Burch,  D.  E.,  J.  Opt.  Soc.  Am.  58,  1383  (1968).  Absorption 
of  Infrared  Radiant  Energy  by  CZS2  H2O.  Ill  Absorption 
of  H2O  between  0.5  and  36  cm“l  (278|j,  -  2cm). 

33.  Burch,  D.  E.,  Gryvnak,  D.  A.,  Patty,  R.  R.,  and  Bartky, 

C.,  E.,  J.  Opt.  Soc.  Am.  59,  267  (1969).  Absorption 

of  Infrared  Radiant  Energ^^y  CO2  and  H2O.  IV.  Shapes 
of  Collision  -  Broadened  CO2  Lines. 

34.  Goulard,  R.,  Editor,  NASA  TMX-53711  (Oct.  1967).  Specialist 
Conference  on  Molecular  Radiation  and  Its  Application  to 
Diagnostic  Techniques. 

35.  Ferrlso,  C.  C.,  J.  Chem.  Phys.  37_,  1955  (1962).  High 
Temperature  Spectral  Absorption  of  the  4.3iJ.  CO2  Band. 

36.  Ferrlso,  C.  C.,  and  Ludwig,  C.  B.,  JQSRT  4,  215  (1964). 
Spectral  Emlsslvltles  and  Integrated  Intensities  of  the 
2.7|jl  H2O  Band  Between  530“  and  2200“K. 

37.  .  Ferrlso,  C.  C.,  and  Ludwig,  C.  B.,  Appl.  Opt.  3,  1435, 

(1964) .  High  Temperature  Spectral  Emlsslvltles  of  H2O, 

CO2  Mixtures  In  the  2.7|i  Region. 

38.  Ferrlso,  C.  C. ,  and  Ludwig,  C.  B.,  J.  Chem.  Phys.  1668, 
(1964)  .  Spectral  Emlsslvltles  and  Integrated  Intensities 
of  the  1.87,  1.38,  1.14|j.  H2O  Bands  Between  1000  and  2200° K. 

39.  Ferrlso,  C.  C.,  Ludwig,  C.  B.,  Thomson,  J.  A.  L. ,  GD/C  - 
DBE  -  65  -  028  (Dec.  1965).  JQSRT  6,  241-273  (1966). 
Empirical  Infrared  Absorption  CoeffTclents  of  H2O  from 
300“K  to  3000“K. 

40.  Ferrlso,  C.  C.,  and  Ludwig,  C.  B.,  Appl.  Opt.  4,  47 
(1965) .  An  Infrared  Band  Ratio  Technique  for  Temperature 
Determinations  of  Hot  Gases. 

41.  Ludwig,  C.  B.,  Ferrlso,  C.  C.,  and  Abeyta,  C.  N.,  JQSRT 
£,  281  (1965).  Emlsslvlty  of  6.3|j,  Band,  H2O  In  Flames. 

42.  Ferrlso,  C.  C.,  Ludwig,  C.  B.,  Thomson,  J.  A.  L.,  Malkmus, 

.  W.,  GD/C  -  DBE  -  66  -  001  (Jan.  1966).  Study  on  Exhaust 

Plume  Radiation  Predictions. 


54 


AEDC-TR-71-78 


43.  Ferriso,  C.  C.,  Ludwig,  C.  B.,  Thomson,  J.  A.  L.,  Malkmus, 

W.,  GD/C  -  DBE  -  66  001a  (Feb.  1966).  Study  on  Exhaust 

Plume  Radiation  Predictions. 

44.  Ludwig,  C.  B.,  Thomson,  J.  A.  L.,  Malkmus,  W.,  GDC  -  DBE 

66- 017  (Dec.  1966),  Final  Report  Contract  NAS  8-11363. 

Study  on  Exhaust  Plume  Radiation  Predictions. 

45.  Ludwig,  C.  B.,  Thomson,  J.  A.  L.,  Malkmus,  W. ,  GDC  -  DBE 

67- 021  (Nov.  1967).  Study  on  Exhaust  Plume  Radiation 
Predictions . 

46.  Ludwig,  C.  B.,  Ferrlso,  C.  C.,  and  Action,  L.,  J.  Opt. 

Soc.  Am.  56,  1685  (1966).  High  Temperature  Spectral 
Emisslvities  and  Total  Intensities  of  the  15-|jl  Band 
System  of  CO2. 

47.  Ludwig,  C.  B.,  and  Ferriso,  C.  C.  JQSRT  7,  7-26,  (1967). 
Prediction  of  Total  Emisslvity  of  Nitrogen  Broadened  and 
Self  Broadened  Hot  Water  Vapor. 

48.  Malkmus,  W.,  and  Thomson,  A.  J.  JQSRT,  2,  17  (1962). 

Infrared  Emisslvity  of  Diatomic  Gases  for  the  Anharmonic 
Vibration  Rotator  Model. 

49.  Malkmus,  W.  J. ,  Opt.  Soc.  Am.  53,  951  (1963).  Infrared 
Emisslvity  of  Carbon  Dioxide  (¥73|jl  Band)  . 

50.  Malkmus,  W.  J.,  Opt.  Soc.  Am.  54,  751  (1964).  Infrared 
Emisslvity  of  Carbon  Dioxide  (2 .7^,  Band)  . 

51.  Malkmus,  W.  J.,  Opt.  Soc.  Am.  57,  323,  (1967).  Random 
Lorentz  Band  Model  with  Exponential  -  Tailed  S~1  line 
Intensity  Distribution  Function. 

52.  Herget,  W.  F.,  Molecular  Radiation  and  Its  Application 
to  Diagnostic  Techniques,  Ed.  by  R.  Goulard  (NASA 
TMX-53711) ,  (1967),  p.  359-379,  Temperature  and  Concen¬ 
tration  Measurements  in  Model  Exhaust  Plumes  Using 
Inversion  Techniques. 

53.  Herget,  W.  F.,  Muirhead,  J.  S.,  and  Golden,  S.  A.,  Report 
No.  R6916,  Rocket dyne  27, (Jan.  1967).  Band  Model  Parameters 
of  H2O . 

54.  Herget,  W.  F.,  and  Muirhead,  J.  S.,  NASA  CR98160; 

Rocketdyne  R7356,  January  1968.  Infrared  Spectral 
Absorption  Coefficient  Data  for  H2O. 

55 .  Herget ,  W .  F . ,  and  Muirhead ,  J .  S . ,  J .  Opt .  Soc .  Am . 

60,  180  (1970).  Infrared  Spectral  Absorption 'Coefficients 
Tor  Water  Vapor. 


55 


AEDC-TR-71-78 


56.  Simmons,  F.  S.,  Appl.  Opt.  ^,1801  (1966).  Band  Models 
for  Nonlsothermal  Radiating  Gases. 

57.  Simmons,  Frederick  S.,  JQSRT,  7,  111  (1967).  Radiances 
and  Equivalent  Widths  of  LorenT'z  Lines  for  Non-isothermal 
Paths . 

58.  Simmons,  F.  S.,'Appl.  Opt.  6,  1423  (1967).  Further  Notes 
on  Band  Models  for  Non-isotliermal  Radiating  Gases. 

59.  Simmons,  F.  S.,  Yamada,  H.  Y. ,  and  Arnold,  C.  B.,  NASA 
CR-72491;  WRL  8962-18-F  (1969).  Infrared  and  Optics 
Laboratory,  Willow  Run  Laboratories,  The  TJniv.  of  Michigan. 
Measurement  of  Temperature  Profiles  in  Hot  Gases  by 
Emission-Absorption  Spectroscopy. 

60.  Tourin,  R.  H.,  J.  Opt.  Soc.  Am.  51,  175  (1961).  Measure¬ 
ments  of  Infrared  Spectral  EmissTvities  of  Hot  Carbon 
Dioxide  in  the  4.3^1  Region. 

61.  Tourin,  R.  H.,  J.  Opt.  Soc.  Am.  799  (1961),  Spectral 

Emisslvlties  of  Hot  CO2-H2O  Mixtures  in  the  2.7iJ,  Region. 

62.  Babrov,  H.  J.,  J.  Opt.  Soc.  Am.  171  (1961).  Instru¬ 

mental  Effects  in  Infrared  Gas  Spectra  and  Spectroscopic 
Temperature  Measurements « 

63.  Ryan,  L.  H.|  Babrov,  H.  J.,  and  Toxurin,  R.  H.,  Warner 
and  Swasey  CO  Report  (1962).  Infrared  Spectra  and 
Temperatures  of  Plasma  Jets. 

64.  Babrov,  H.  J.,  J.  Opt.  Soc.  Am.  53,  945  (1963).  Experi¬ 

mental  and  Theoretical  Spectral  )[Fsorptance  of  HCl  at 
Various  Temperatures. 

65.  Babrov,  H.  J.,  and  Tourin,  R.  H. ,  JQSRT  3,  15-27,  (1963). 
Methods  for  Predicting  Infrared  Radiance  of  Flames  by 
Extrapolation  of  Laboratory  Measurements. 

66.,  Tourin,  R.  H.,  and  Krakow,  B.,  Appl.  Opt.  4,  237,  (1965). 
Applicability  of  Infrared  Emission  and  Absorption 
Spectra  to  Determination  of  Hot  Gas  Temperature  Profiles. 

67.,  Krakow,  B.',  Tourin,  R.  H.,  Penzias,  G.  J.,  NASA  CR-66 112, 
Warner  and  Swasey  CO  Report  TR-39.  Possible  Techniques 
for  Optical  Measurement  of  Temperature  and  Concentration 
Profiles  in  a  Supersonic  Ramjet. 

68.  Krakow,  B.,  Appl.  Opt.  5,  201  (1966).  Spectroscopic 
Temperature  Profile  Measurements  in  Inhomogeneous  Hot  Gases. 

69.  Krakow,  B.,  Babrov,  H.  J.,  Maclay,  G.  J.,  and  Shabott, 

A.  L.,  Appl.' Opt.  5,  1791  (1966).  Use  of  the  Curtis- 
Godson  Approximation  in  Calculations  of  Radiant  Heating 
by  Inhomogeous  Hot  Gases. 


56 


AEDC-TR-71-78 


70.  Babrov,  H.  J.,  Shabott,  A.  L. ,  and  Krakow,  B.,  Final 
Report,  Part  3  (TR-41) .  Contract  NAS  8-11315  (Mar.  1966). 
Applications  of  Spectroscopic  Absorption  Data  and 
Inhomogeneous  Transmittance  Combination  Methods. 

71.  Tourin,  R.  H.,  "Spectroscopic  Gas  Temperature  Measurement; 
Pyrometry  of  Hot  Gases  and  Plasmas",  Elsevier  Publishing 
Company,  N.  Y. ,  (1966). 

72.  Krakow,  B.,  and  Babrov,  H.  J.,  Molecular  Spectral 
Parameters  for  Correlation  of  Spectral  Radiance  with  Gas 
Temperature  Profiles. 

73.  Babrov,  H.  J.,  and  Casden,  F.  J.  Opt.  Soc.  Am.  58,  179 

(1968).  Strengths  of  Forty-two  Lines  in  the  V3  Bands 

of  Water  Vapor. 

74.  Babrov,  H.  J.,  Healy,  A.  R.,  J.  Opt.  Soc.  Am.  779 

(1969) .  Strengths  of  Twenty  Lines  in  the  V3  Band  of 
Water  Vapor. 

75.  Krakow,  B.,  and  Healy,  A.  R.,  J.  Opt.  Soc.  Am.  1490, 

(1969) .  Strengths  of  31  Water  Vapor  Lines  Between  1617 
and  1429  cm-1. 

76.  Oppenheim,  U.  P. ,  and  Ben-Aryeh,  JQSRT  4,  559-570, 

(1964).  A  General  Method  for  the  Use  oT  Band  Models, 
with  Application  to  Infrared  Atmospheric  Absorption. 

77.  Oppenheim,  U.  P.,  and  Ben-Aryeh,  J.  Opt.  Soc.  Am.  53, 

344-350,  (1963)  .  Statistical  Model  Applied  to  the'TJegion 
of  the  V3  Fundamental  of  CO2  a-t  1200°K. 

78.  Y.  Ben-Aryeh,  Appl.  Opt.,  6,  1049-1055  (1967).  Spectral 
Emissivity  Calculations  by  the  Statistical  Model  Applied 
to  the  4.3p,  Bands  of  CO2  a-t  High  Temperatures. 

79.  Ellis,  R.  E.,  and  Schurin,  B. ,  Appl.  Opt.  8,  2265  (1969). 
Integrated  Intensity  Measurements  of  CO2  Bands  in  the 
4.82p.  and  5.17|j.  Regions. 

80.  Kaplan,  L.  D.,  and  Eggers,  D.  F.,  J.  Chem.  Phys.  876 

(1956)  .  Band  Strength  of  15^1  CO2  Band  Using  Curve  of 
Growth  (nitrogen  broadening) . 

81.  Varanasi,  P.,  and  Lauer,  J.  L.,  JQSRT  6,  127-130,  (1966). 

A  Re  determination  of  the  Integrated  Intensity  of  the  15|i. 
Bands  of  Carbon  Dioxide . 

82.  Wolk,  M.,  JQSRT  7,  1-5,  (1967).  The  Strength  of  the  Pressure 
Broadened  CO2  Bands  at  15|x  by  Digital  Integration  of  Spectra. 


57 


AEDC-TR-71-78 


83.  Gray,  L.  D.,  JQSRT  7,  143  (1967).  Calculations  of  Carbon 
Dioxide  Transmission.  Part  1,  the  9.4p.  and  10.4|x  Bands. 

84.  Harward,  D.  N.,  and  Patty,  R.  R. ,  J.  Opt.  Soc.  Am.  58 , 

188  (1968) .  Low  Resolution  Determination  of  the  Strength 
of  the  667-cm“l  CO2  Band. 

85.  Gates,  D.  M. ,  Calfee,  R.  F.,  Hanson,  D.  W.,  and  Benedict, 

W.  S.,  National  Bureau  of  Standards  Monograph  71,  (1964) 

Line  Parameters  and  Computed  Spectra  for  Water  Vapor 
Bands  at  2.7|jl. 

86.  Benedict,  W.  S.,  and  Calfee,  R.  F.,  ESSA  Professional 
Paper  2,  U.  S.  Department  of  Commerce,  (1967).  Line 
Parameters  for  the  1.9  and  6.3  Micron  Water  Vapor  Bands. 

87.  Benedict,  W.  S.,  and  Kaplan,  L.  D.,  JQSRT  4,  453  (1964). 
Calculation  of  Line  Widths  in  H2O-H2O  and  E2O-O2  Collisions. 

88.  Goldman,  A.,  and  Oppenhelm,  U.  P. ,  J.  Opt.  Soc.  Am.  55, 

794  (1965) .  Emissivity  of  Water  Vapor  at  1200®K.  in 
the  1.9  -  and  2.7-|i  Regions. 

89.  Hlrshfeld,  M.  A.,  Jaffe,  J.  H. ,  and  Ross,  G. ,  JQSRT  6, 
311-315  (1966).  New  Determination  of  the  Strength  of  ■ 
the  U3  Band  of  Water  Vapor  From  Dispersion  Measurements. 

90.  Aryeh,  Y.  Ben,  JQSRT  7,  211  (1967).  Line  Widths  and 
Intensities  in  the  Wings  of  the  Up  Water  Vapor  Band  at 
400®K  and  540®K. 

91.  Izatt,  J.  R.,  Sakai,  H.,  and  Benedict,  W.  S.,  J.  Opt. 

Soc.  Am.,  59,  19,  (1969).  Positions ,' Intensities ,  and 
Widths  of  Walter  Vapor  Lines  Between  475  and  692  cm"l. 

92.  Fraley,  P.  E.,  Rao,  K.  N.,  and  Jones,  L.  H.,  J.  Mol. 

Spectry  29,  312  (1969).  High  Resolution  Infrared  Spectra 
of  Water"^apor  vi  and  V3  Bands  of  H2^®0. 

93.  Fraley,  P.  £.,  and  Rao,  K.  N. ,  J.  Mol.  Spectry  29,  348, 
(1969) .  High  Resolutign  Infrared  Spectra  of  WaTer  Vapor 
Vj  and  V3  Bands  of  H2^®0. 

94.  Stull,  R.  V.,  and  Plass ,  G.  N.,  J.  Opt.  Soc.  Am.,  50, 

1279  ^960).  Band  Model:  Emissivity  of  HCl. 

95.  Gryvnak,  D.  A.,  and  Shaw,  J.  H.,  J.  Opt.  Soc.  Am.,  52, 

539  (1962).  Study  of  the  Total  Absorption  Near  4.7)i 
by  Two  Samples  of  CO  as  Their  Total  Pressures  and  CO 
Concentrations  were  Independently  Varied. 


58 


AEDC-TR-71-78 


96.  Abu-Romia,  M.  M.,  and  Tien,  C.  L. ,  JQSRT  6,  143-167, 

(1966) .  Measurements  and  Correlations  of  infrared 
Radiation  of  Carbon  Monoxide  at  Elevated  Temperatures 
(300,  600,  900,  1200,  1500®K) (^  to  3  atm). 

97.  Gray,  L.  D.,  and  McClatchey,  R.  A.,  Appl.  Opt.  4,  1624 
(1965).  Calculations  of  Atmospheric  Radiation  From 
4.2^  to  5^.  C02f  NoO,  CO,  Random  Elsasser  Band  Model 
with  Curtis-Godson  Approx,  for  Inhomogeneous  Paths. 

98.  Gray,  L.  D.,  Appl.  Opt.  4,  1494  (1965).  Spectral 
Absorption  of  the  4.6^  Bands  of  N2O. 

99.  Jacobs,  T.  A.,  and  Giedt,  R.  R.,  JQSRT  3,  289-291  (1963). 
The  Applicability  of  Beer's  Law  for  a  Broad  Spectral 
Bandpass . 

100.  Golden,  S.  A.,  The  Doppler  Analog  of  the  Elsasser  Band 
Model,  Contract  AF  33(615)  3093. 

101.  Armstrong,  B.  H. ,  JQSRT  7,  61-88,  (1967).  Spectrum  Line 
Profiles:  The  Voigt  Function,  A  Review  of  Math  Properties, 

.  Recent  Literature,  Computational  Methods. 

102.  Whiting,  E.  E.,  Arnold,  J.  O. ,  Lyle,  G.  C.,  NASA  TND  - 
5088,  (March  1969).  A  Computer  Program  for  a  Line-by- 
Line  Calculation  of  Spectra  from  Diatomic  Modecules  and 
Atoms  Assuming  a  Voigt  Line  Profile. 

103.  Kyle,  T.  G.,  J.  Opt.  Soc.  Am.  58,  192  (1968).  Net  Interval 
for  Calculating  Absorption  SpecTra. 

104.  Watson,  R.,  Planet,  W.  G.,  and  Pitts,  C.  C.,  Appl.  Opt.  7 
1941  (1968) .  On  Temperature  Determinations  from  Non- 
resolved  Spectra.  Discussion  of  Instrumental  Slit 
Function  Effects  in  the  Spectrograph-film-densitometer 
Combination. 

105.  Goldman,  A.,  and  Kyle,  T.  G.. ,  Appl.  Opt.  7,  1167  (1968). 

A  Comparison  Between  Statistical  Model  anH*  Line  by  Line 
Calculation  with  Application  to  the  9.6|x  Ozone  and  the 
2.7iJ,  Water  Vapor  Bands. 

106.  Rusk,  J.  R.,  J.  Chem.  Phys.  42,  493  (1965),  Line  Breadth 
Study  of  the  1.64  mm  Absorption  in  Water  Vapor. 

107.  Silverman,  S.,  J.  Opt.  Soc.  Am.  39,  275  (1949). 


59 


60 


AEDC-TR-71-78 


Fig.  4  Geometry  of  Lines  of  Sight  for  Radially 
Symmetric  Flame 
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PI.  -  Plexiglass  Enclosures  S.M.  -  Spherical  Mirror  H.  -  Exhaust  Hood 

P.E.  -  Purge  Gas  Entrance  (Ar)  P.M.  -  Plane  Mirror  M.  -  Monochromator 

G.  -  Globar  (Continuum  Source)  B.  -  Burner  (1-1/16"  diam.)  -  Absorption  Optical 

C.  -  Chopper  (300  cps)  Sh.  -  Shutter  (Employed  for  Path 

S-,  -  Entrance  Slit  Emission  Only)  .  Emission  Optical 

^  Path 

a  -  20  deg 


Fig.  5  Schematic  of  the  Optical  System 
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Radiance  (watts/ster-cm  -ij.) 


Fig.  6  Measured  Spectral  Radiance  and  Transmissivity  Profiles  for  a 
Hydrogen-Air-Carbon  Dioxide  {4%)  Flame 
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Fig.  7  Temperature  and  Water  Vapor  Partial  Pressure  Profiles  for 
Hydrogen-Air-Carbon  Dioxide  (4%)  Flame  (Neglecting 
Foreign  Gas  Broadening) 
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TABLE  I 

STATION  0.  HEIGHT  ABOVE  BURNER  BASE  1/2" 
Minus  side  of  flame 


Zone 

«i 

« 

Th 

Ph 

1 

0.0088 

0.9946 

1329.768 

0.160144 

2 

0.0235 

0.9908 

1598.593 

0.1912517 

3 

0.0439 

0.9884 

1830.329 

0.2416925 

4 

0.0683 

0.984* 

1766.159 

0.3725197 

5 

0.0964 

0.9816 

2042.532 

0.3851425 

6 

0.113 

0.9822 

2415.941 

0.2547318 

7 

0.126 

0.9802 

1973.946 

0.4243071 

8 

0.1334 

0.9788 

1976.341 

0.3987573 

9 

0.134 

0.9782 

1868.96 

0.3734008 

10 

0.134 

0.9761 

1425.101 

0.7195553 

Plus  side  of  flame 

Zone 

«i 

*^1 

Ph 

1 

0.0072 

0.9929 

1199.679 

0.1040097 

2 

0.0202 

0.9929 

1633.923 

0.1072218 

3 

0.0439 

0.9891 

1765.506 

0.1960801 

4 

0.0675 

0.9864 

1881.406 

0.2533221 

5 

0.0952 

0.9826 

1914.978 

0.3425443 

6 

0.1098 

0.9812 

1961.568 

0.3467113 

7 

0. 1204 

0.9795 

1939.124 

0.3689625 

8 

0.1286 

0.9787 

1968.003 

0.3739443 

9 

0.134 

0.9778 

1950.502 

0.3966639 

10 

0.134 

0.9761 

1743.325 

0.4711858 
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TABLE  II 

STATION  1.  HEIGHT  ABOVE  BURNER  BASE  3/4" 
Minus  side  of  flame 


Zone 

H. 

1 

Ph 

1 

0.004 

0.9959 

1190.162 

0.06572864 

2 

0.0138 

0.9931 

1522.641 

0.1345853 

3 

0.0333 

0.9917 

2039.2 

0.1599885 

4 

0.0577 

0. 9886 

1914.87 

0.2729758 

5 

0.0817 

0.9823 

1704.095 

0.4321557 

6 

0.109 

0.9819 

2399.626 

0.2897842 

7 

0.12 

0. 9805 

1918.579 

0.3646834 

8 

0.137 

0.9788 

2116.154 

0.4257887 

9 

0.14 

0. 9788 

1870.317 

0.3706299 

10 

0.14 

0.9778 

2008.301 

0.3253451 

11 

0.138 

0.9784 

2556.197 

0.1313624 

Plus  side  of  flame 

Zone 

«1 

^1 

Ph 

1 

0.004 

0.9948 

1134.789 

0.07272909 

2 

0.0162 

0.9917 

1434.445 

0.1245156 

3 

0.0389 

0.9889 

1710.057 

0.1905522 

4 

0.0667 

0.9826 

1700.995 

0.3211193 

5 

0.0928 

0.9812 

1898.645 

0. 3441612 

6 

0.1098 

0. 9809 

2045.449 

0.3298195 

7 

0. 1254 

0.9802 

2132.306 

0.3394778 

8 

0.1344 

0.9798 

2172.490 

0.3376801 

9 

0. 1376 

0.9784 

2005.149 

0.3855848 

10 

0.1384 

0.9784 

2027.117 

0.3704518 
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TABLE  III 

STATION  2.  HEIGHT  ABOVE  BURNER  BASE  1" 
Minus  side  of  flame 


Zone 

«1 

”^1 

Ph 

i  ' 

0.0040 

0.9932 

1078.085 

0.1008756 

2 

0.0121 

0.9889 

1308.032 

0.1804165 

3 

0.0292 

0.9864 

1667.921 

0.2097130 

4 

0.0517 

0.9846 

1883.424 

0.2570956 

5 

0.0777 

0.9828 

2048.448 

0.3065125 

6 

0.1017 

0.9806 

2041.374 

0.3651089 

7 

0.1266 

0.9799 

2386.426 

0.3304918 

8 

0.1353 

0.9792 

2126.313 

0.3350262 

9 

0. 1388 

0.9774 

1874.469 

0.4126021 

10 

0.1416 

0.9767 

1957.927 

0.3834246 

11  . 

0.1418 

0.9765 

1849.66 

.  0.3792502 

Plus  side  of  flame 


Zone 

Ph 

1  ■ 

0.0113 

0.995 

1436.612 

0.08800917 

2 

0.0276 

0.9918 

1623.908 

0.1455148 

3 

0.0483 

0.9903 

1879.595 

0.1799156 

4 

0.0737 

0.9889 

2110.551 

0.2159047 

5 

0.1017 

0.9853 

2059.729 

0.3108270 

6 

0.1196 

0.9849 

2229.696 

0.2931871 

7 

0.1298 

0.9824 

2047.16 

0.3532553 

8  ' 

0.1364 

0.9799 

1907.419 

0.4071964 

9 

0.1392 

0.9788 

1865.414 

0.4166042 

10 

0. 1408 

0.9785 

1869.995 

0.4215534 
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TABLE  IV 

STATION  3.  HEIGHT  ABOVE  BURNER  BASE  r/4" 
Minus  side  of  flame 


Zone 

Ph 

1 

0.008 

0.9942 

1281.014 

0.09960471 

2 

0.017 

0.9893 

1348.209 

0.1967253 

3 

0.03 

0.986 

1536.079 

0.2175365 

4 

0.0521 

0.9847 

1933.05 

0.2439525 

5 

0.0781 

0.984 

2212.649 

0.2597838 

6 

0.1046 

0.9821 

2157.58 

0.3459801 

7 

0.1242 

0.9808 

2227.028 

0.3393879 

8 

0.1364 

0.9788 

2005.548 

0.4149947 

9 

0. 1404 

0.9785 

2130.546 

0.3201958 

10 

0.142 

0.9778 

1919.647 

0.3823121 

11 

0.1426 

0.9762 

1535.3^ 

0.6208833 

Plus  side  of  flame 

Zone 

*^1 

1 

0.0088 

0.9883 

1129,815 

0.1636166 

2 

0.0243 

0.988 

1487.810 

0.1630804 

3 

0.0475 

0.984 

1638.096 

0.2439516 

4 

0.0765 

0.9834 

1957.051 

0.2689019 

5 

0. 1005 

0.9821 

2084.923 

0.2968581 

6 

0.1148 

0.9817 

2177.459 

0.2881885 

7 

0.1262 

0.9804 

2127.871 

0.3184034 

8 

0.1348 

0.9798 

2153.448 

0.3254421 

9 

0.1384 

0.9765 

1826.465 

0.4465935 

10 

0.1416 

0.9762 

1873.864 

0.4548999 
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